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出版说明 

■HHBIBIBMMMMQaMMBiaHaiMMIRRMnMMMHMBinMHIMmHnR 


由香港科技大学社会科学部吴晓刚教授主编的“格致方 
法 • 定量研究系列”丛书，精选了世界著名的 SAGE 岀版社 
定量社会科学研究丛书，翻译成中文，起初集结成八册，于 
2011年岀版。这套丛书自出版以来，受到广大读者特别是年 
轻一代社会科学工作者的热烈欢迎。为了给广大读者提供 
更多的方便和选择，该丛书经过修订和校正，于2012年以单 
行本的形式再次出版发行，共37本。我们衷心感谢广大读 
者的支持和建议。 

随着与 SAGE 岀版社合作的进一步深化，我们又从丛书 
中精选了三十多个品种，译成中文，以繪读者。丛书新增品 
种涵盖了更多的定量研究方法。我们希望本丛书单行本的 
继续岀版能为推动国内社会科学定量研究的教学和研究作 
岀一点贡献。 




2003 年，我赴港工作，在香港科技大学社会科学部教授 
研究生的两门核心定量方法课程。香港科技大学社会科学 
部自创建以来，非常重视社会科学研究方法论的训练。我开 
设的第一门课“社会科学里的统计学 ” （Statistics for Social 
Science ) 为所有研究型硕士生和博士生的必修课，而第二门 
课“社会科学中的定量分析”为博士生的必修课（事实上，大 
部分硕士生在修完第一门课后都会继续选修第二门课）。我 
在讲授这两门课的时候，根据社会科学研究生的数理基础比 
较薄弱的特点，尽量避免复杂的数学公式推导，而用具体的 
例子，结合语言和图形，帮助学生理解统计的基本概念和模 
型。课程的重点放在如何应用定量分析模型研究社会实际 
问题上，即社会研究者主要为定量统计方法的“消费者”而非 
“生产者”。作为“消费者”，学完这些课程后，我们一方面能 
够读懂、欣赏和评价别人在同行评议的刊物上发表的定量研 
究的 文章; 另一方面，也能在自己的研究中运用这些成熟的 
方法论技术。 

上述两门课的内容，尽管在线性回归模型的内容上有少 
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量重复，但各有侧重。“社会科学里的统计学”从介绍最基本 
的社会研究方法论和统计学原理开始，到多元线性回归模型 
结束，内容涵盖了描述性统计的基本方法、统计推论的原理、 
假设检验、列联表分析、方差和协方差分析、简单线性回归模 
型、多元线性回归模型，以及线性回归模型的假设和模型诊 
断。“社会科学中的定量分析”则介绍在经典线性回归模型 
的假设不成立的情况下的一些模型和方法，将重点放在因变 
量为定类数据的分析模型上，包括两分类的 logistic 回归模 
型、多分类 logistic 回归模型、定序 logistic 回归模型、条件 lo ¬ 
gistic 回归模型、多维列联表的对数线性和对数乘积模型、有 
关删节数据的模型、纵贯数据的分析模型，包括追踪研究和 
事件史的分析方法。这些模型在社会科学研究中有着更加 
广泛的应用。 

修读过这些课程的香港科技大学的研究生，一直鼓励和 
支持我将两门课的讲稿结集岀版，并帮助我将原来的英文课 
程讲稿译成了中文。但是，由于种种原因，这两本书拖了多 
年还没有完成。世界著名的出版社 SAGE 的“定量社会科学 
研究”丛书闻名遐迩，每本书都写得通俗易懂，与我的教学理 
念是相通的。当格致出版社向我提出从这套丛书中精选一 
批翻译，以飨中文读者时，我非常支持这个想法，因为这从某 
种程度上弥补了我的教科书未能出版的遗憾。 

翻译是一件吃力不讨好的事。不但要有对中英文两种 
语言的精准把握能力，还要有对实质内容有较深的理解能 
力，而这套丛书涵盖的又恰恰是社会科学中技术性非常强的 
内容，只有语言能力是远远不能胜任的。在短短的一年时间 
里，我们组织了来自中国内地及香港、台湾地区的二十几位 



研究生参与了这项工程，他们当时大部分是香港科技大学的 
硕士和博士研究生，受过严格的社会科学统计方法的训练， 
也有来自美国等地对定量研究感兴趣的博士研究生。他们 
是香港科技大学社会科学部博士研究生蒋勤、李骏、盛智明、 
叶华、张卓妮、郑冰岛，硕士研究生贺光烨、李兰、林毓玲、肖 
东亮、辛济云、於嘉、余珊珊，应用社会经济研究中心研究员 
李 俊秀; 香港大学教育学院博士研究生洪 岩璧; 北京大学社 
会学系博士研究生李丁、赵亮员；中国人民大学人口学系讲 
师巫 锡炜； 中国台湾“中央”研究院社会学所助理研究员林宗 
弘•，南京师范大学心理学系副教授陈陈；美国北卡罗来纳大 
学教堂山分校社会学系博士候选人姜 念涛; 美国加州大学洛 
杉矶分校社会学系博士研究生 宋曦; 哈佛大学社会学系博士 
研究生郭茂灿和周韵。 

参与这项工作的许多译者目前都已经毕业，大多成为中 
国内地以及香港、台湾等地区高校和研究机构定量社会科学 
方法教学和研究的骨干。不少译者反映，翻译工作本身也是 
他们学习相关定量方法的有效途径。鉴于此，当格致出版社 
和 SAGE 岀版社决定在“格致方法 • 定量研究系列”丛书中 
推出另外一批新品种时，香港科技大学社会科学部的研究生 
仍然是主要力量。特别值得一提的是，香港科技大学应用社 
会经济研究中心与上海大学社会学院自2012年夏季开始， 
在上海(夏季）和广州南沙(冬季)联合举办《应用社会科学研 
究方法研修班》，至今已经成功举办三届。研修课程设计体 
现“化整为零、循序渐进、中文教学、学以致用”的方针，吸引 
了一大批有志于从事定量社会科学研究的博士生和青年学 
者。他们中的不少人也参与了翻译和校对的工作。他们在 
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繁忙的学习和研究之余，历经近两年的时间，完成了三十多 
本新书的翻译任务，使得“格致方法 • 定量研究系列”丛书更 
加丰富和完善。他们是：东南大学社会学系副教授洪岩璧， 
香港科技大学社会科学部博士研究生贺光烨、李忠路、王佳、 
王彦蓉、许多多，硕士研究生范新光、缪佳、武玲蔚、臧晓露、 
曾东林，原硕士研究生李兰，密歇根大学社会学系博士研究 
生王骁，纽约大学社会学系博士研究生温芳琪，牛津大学社 
会学系研究生周穆之，上海大学社会学院博士研究生陈 
伟等。 

陈伟、范新光、贺光烨、洪岩璧、李忠路、缪佳、王佳、武玲 
蔚、许多多、曾东林、周穆之，以及香港科技大学社会科学部 
硕士研究生陈佳莹，上海大学社会学院硕士研究生梁海祥还 
协助主编做了大量的审校工作。格致出版社编辑高璇不遗 
余力地推动本丛书的继续出版，并且在这个过程中表现出极 
大的耐心和高度的专业 精神。 对他们付出的劳动，我在此致 
以诚挚的谢意。当然，每本书因本身内容和译者的行文风格 
有所差异，校对未免挂一漏万，术语的标准译法方面还有很 
大的改进空间。我们欢迎广大读者提出建设性的批评和建 
议，以便再版时修订。 

我们希望本丛书的持续出版，能为进一步提升国内社会 
科学定量教学和研究水平作岀一点贡献。 


吴晓刚 

于香港九龙清水湾 



标注说明 


黑体字母用来表示矩阵和向量，例如 B 表明 B 是一个 
矩阵。矩阵和向量的维度用 ( rXc ) 来表示^代表行的数量 ， c 
代表列的数量。用带有下标的小写字母来表示矩阵的元素， 
例如 B 矩阵中的第 ( ij ) 个元素为 b ; j ， x 向量的第 i 个元素为: Tj 。 

表示转置矩阵，例如以是 B 的转置矩阵。带有上角标 “一1” 
的矩阵用来表示逆矩阵，如 B - 1 为 B 的逆矩阵。 “ COV ” 表示 
协方差，如果括号里面为两个变量，则 COV ( x ,, : rj 表示两 
个变量之间的协方差，如果括号里是向量集合，比如 （nX 1) 
维的 x 集合，那么 COV ( x ) 就是 （nX n ) 维的协方差矩阵，其 
中第 ( ij ) 个元素 ( i 关 j ) 表示 A 和: Tj 之间的协方差，第 （ ii ) 个 
元素表示： Ti 的方差。同样, “ COR ” 表示相关， COR (:^， &) 
表示两个变量 x , 和&之间的 相关; COR ( x ) 则为 （ nXii ) 维 
的相关矩阵，其中第 ( ij ) 个元素 ( i 关 j ) 表示^和: Tj 之间的相 
关，第 （ ii ) 个元素为1。 “ E ” 表示数学期望，如果^是随机变 
量，则 ECn ) 为变量^的期 望值; 如果 x 是个向量集，则 E ( x ) 
表示随机变量 A 的期望值中的第 i 个元素。 

本书中的图标、方程、例子、表格等都按照章节的顺序来 
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命名，如表 2.3 表示第2章中的第三 张表; 例 3.2 表示第3章 
中的第二个例子。需要注意的是，本书的例子可能是逐步延 
伸和发展的，如果在第3章中看到有多处例3.2,读者应该注 
意，它们指的是同一个例子。 



《协方差结构模型: LISREL 导论》直接建立在作者《验证 
性因子分析 》 (Conf irmatory Factor Analysis ) 的基础上。建 

议不熟悉验证性因子分析模型的读者先行阅读《验证性因子 
分析》。掌握验证性因子分析模型和矩阵的相关知识有助于 
理解协方差结构模型。 

在导论中，作者简单地介绍了测量模型、结构方程模型 
及整合了两者的协方差结构模型。在第2章中，作者在验证 
性因子模型的基础上讨论了测量模型。对于那些已经掌握 
了《验证性因子分析》的读者而言，本章的内容相当明了，应 
该不会遇到什么困难，只需稍加留意方程和参数的标注 
即可。 

第3章介绍了结构方程模型。阅读本章时，读者需要掌 
握一些计量经济学中关于结构方程模型的知识。为了教学 
的方便，本章在介绍结构方程模型时，假定所有的变量都不 
含测量误差。作者会在其后的章节里放松这一假定，并将测 
量模型和结构方程模型结合起来。作者之所以分别介绍两 
种传统的定量研究方法（即心理学中的测量模型和计量经济 
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学中的结构方程模型），是因为多数读者或许熟悉其中的一 
种方法，少数读者或许对这两种方法都比较熟悉。 

从某种程度上来说，本书的第3章与作者的《验证性因 
子分析》是相互对应的。本章中，作者讨论了各种类型的结 
构方程模型、每种模型的识别和估计方法，以及关于模型拟 
合度的评价等。 

第4章介绍了协方差结构方程模型。对于那些已经熟 
悉《验证性因子分析》和本书第3章的读者来说，本章仅仅是 
一篇综述。相对于验证性因子分析模型，协方差结构分析考 
虑到了潜变量之间的结构 关系； 而相对于结构方程模型，协 
方差结构分析意识到了社会科学中的大部分概念是无法被 
直接观测到的，因而整合了测量模型。虽然整合了两者的协 
方差结构分析增加了模型识别和估计的复杂性，但这更符合 
实际情况。本书还有助于拓宽读者的视野，并向读者介绍了 
心理学和计量经济学相互结合的丰硕成果。 

约翰 • L . 沙利文 



本书所介绍的统计模型可以被称做协方差结构模型、协 
方差结构分析、线性结构模型、矩阵结构模型、线性结构中的 
潜变量方程组，或常见的 LISREL 模型。本书采用协方差结构 
模型这一更为广泛的概念。协方差结构模型的估计需要复杂 
的统计软件，由约瑞斯科 ( Joreskog ) 和索尔波姆 ( S 6 rbom ) 所 
发展的 LISREL 是目前应用最为广泛的协方差结构分析软 
件。 LISREL 的重要性不仅体现在其软件功能上，而且这一词 
语已经成为代表数据分析的一种模型和方法。本书的副标题 
“ LISREL 导论”也体现了其重要性。需要注意的是，本书并非 
关于如何使用 LISREL 软件的人门读物。 

许多读者可能不熟悉复杂的协方差结构模型，但是他们 
或许已经掌握了这个模型的部分知识。因为协方差结构分 
析是由众所周知且极为重要的两种统计技术组成的。第一 
部分是心理学家经常使用的验证性因子 分析; 第二部分则是 
计量经济中的结构方程模型。读者如果至少熟悉其中的一 
种方法，就会比较容易理解本书的内容。 

作者假定读者已经熟悉了《验证性因子分析》中所介绍 
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的验证性因子分析模型及数学工具。对结构方程模型感兴 
趣的读者也许会发现，本书的讨论对于参数等同限定及误差 
相关的估计会很有帮助。 

全面理解协方差结构模型需要将其运用到实际的数据 
分析中，因此，我们鼓励读者运用本书附录中的数据去练习 
和再现书中的例子。如果能得到相似的结果，就表明读者已 
经很好地理解了本书的内容。 一 般的统计软件（如 SPSS 、 
SAS 、 BMDP ) 通常没有关于协方差结构模型的估计程序，本 
书结尾会简单介绍一些可以应用于协方差结构分析的统计 
软件。 

许多人对本书的岀版给予了慷慨的建议和评价。我想 
在此表达对他们的感谢。他们是保罗 • 艾利森 （Paul 
Allison )、 格雷格 • 邓肯 （Greg Duncan )、 卡伦 • 普格利西 
(Karen Pugliesi ) 、杰伊 • 斯图尔特 （Jay Stewart )、布莱尔 • 
惠顿 （Blair Wheaton ) 、罗纳德 • 舍恩伯格 （Ronald Schoen - 
berg )， 及两位匿名评阅人。卡罗尔 • 赫克曼 （Carol 
Heckman ) 通读了本书和《验证性因子分析》的草稿，她的评 
阅极大地提高了本书的精确性和清晰性。当然，本书中的错 
误与上面提到的诸位无关。 
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协方塞结构 横型: LKRKL 导论 


协方差结构分析模型试图用比较少的潜变量来解释一 
组观测变量之间的关系。正如这个方法的名字所暗示的那 
样，观测变量之间关系的是包含在其协方差矩阵5：中的。协 
方差结构分析假定观测变量之间的关系和特征是由一些潜 
变量决定的，通过测量模型将观测变量和潜变量连接起来， 
然后用结构模型来表示潜变量之间的关系。协方差矩阵 S 
的相关分析就是用以描述这种结构的。 

“协方差结构分析”这个术语是由博克 （ Bock ) 和巴格曼 
( Bargmann ) (1996) 引入的，他们曾用这个术语来描述现在被 
称做验证性因子分析的模型。自那以后，许多学者都为提升 
这个模型的综合性和普遍性作出了巨大贡献。这个模型从 
最初由博克和巴格曼引入的因子分析模型，经由很多中间模 
型的发展而成为极具普遍性的模型。在这个模型里，协方差 
矩阵2：可以被看做任何参数的任何函数形式。关于这些中 
间模型的综述可以参阅本特勒和威克斯的论述 （Bemler & 
Weeks , 1979) 或者是本特勒 （ Bender ，1980) 等学者的论述。 

尽管协方差结构模型在估计和应用上都有了长足的发 
展，但本书重点介绍由约瑞斯科 ( hreskog ， 1973)、约瑞斯科 
和范西洛 （ j 6 reskog van Thillo ，1972) 、基斯林 （ Keesling ， 
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1972) 和威利 ( Wiley , 1973) 等人引人的那些较为简洁但又普 
遍应用的模型。这些简洁模型将观测变量的协方差分解为 
概念上相互独立的两个部分。第一步，类似于心理测量学的 
做法，通过因子分析模型将观测变量和潜变量连接起来。第 
二步，类似于计量经济学的做法，用结构模型来设定潜变量 
之间的因果关系。协方差结构模型整合了心理测量学和计 
量经济学，可以同时设定因子模型和结构模型。这一整合很 
大程度上得益于戈德伯格 （ Goldberger ) 于1970年组织的结 
构方程模型会议 (Goldberger 8^ Duncan ，1973) 及其1971年 
关于结构方程模型程序的文章。 

协方差结构模型的应用需要使用能够使得包含众多变 
量的函数方程最大化的高效的数值方法。1966年，当时还在 
教育服务机构工作的约瑞斯科在这一领域作岀了重大突破。 
一系列的程序被开发出来，最终发展成为著名的且被广泛使 
用于协方差结构模型分析的 LISREL 软件 ( J 6 resk 0 g & van 
Thillo , 1972； Joreskog Sorbom , 1976，1978， 1981)。 
LISREL 为协方差结构模型的普及和应用作岀了巨大的贡 
献，以致协方差结构模型常常被称做 “ LISREL 模型”。 

本书中的协方差结构模型实际上是指那些被整合进 
LISREL 的模型，下文统称“协方差结构模型”。我们之所以 
只关注 LISREL 中的协方差结构模型是基于以下三点 考虑: 
第一， LISREL 中有许多程序可以估计协方差结构模型。这 
很重要，因为如果没有估计软件，我们很难应用这些模型。 
第二，更广泛的协方差结构模型对数学的要求超过了本书的 
范围。第三，这些模型包括了验证性因子分析、二阶因子分 
析，多指标模型、联立方程模型、追踪数据模型，以及方程和 
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变量中含有误差项的结构方程模型等，这些模型应该可以满 
足大部分研究者的需求。 

在本书中，我们首先分别介绍因子分析模型和结构方程 
模型，然后将它们整合到一起再来介绍协方差结构模型。第 
2章主要介绍协方差结构模型中的测量部分，这和验证性因 
子分析模型是相似的。第3章主要介绍协方差结构分析中 
的结构方程模型，它可以被看做只包含结构部分的协方差结 
构分析的一个特例。第4章的介绍同时整合了第2章中的 
因子模型和第3章中的结构模型的协方差结构模型。 

下面，我们首先简要介绍每个模型中的数学结构。 

因子分析模型假定观测变量是由较少的潜变量生成的， 
这些潜变量被称做因子。由于观测变量通常会含有测量误 
差，因此因子模型基本上可以被看做测量模型。例如图 1.1 
中的模型，方框对应的是含有测量误差的观测变量，圆圈对 
应的是潜变量。顶部的两个圆圈中的潜变量分别对应着两 
个观测变量，表示这些观测变量是由潜变量产生的，这两个 
潜变量被称做共同因子。底部的每个圆圈分别对应一个观 
测变量，这些圆圈表示观测变量中不能由共同因子解释的部 
分，因此它们被称做独特因子或者测量误差。尽管可以知道 
观测变量之间的协方差，可是它们已经受到了测量误差的干 
扰。为了移除测量误差的影响，测量模型只估计图 1.1 中顶 
部两个共同因子之间的协方差。 

图 1.1 中潜变量和观测变量之间的关系可以通过下列方 
程来 表示： 

工 1 ~Anfi d\ JCz = A2i6i + 占 2 

工 3 =又32?2 + 谷 3 4 ~ A 42 f 2 H - ^4 [1.1] 
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图 U 协方差结构模型中的测置部分 


方程中的: r 是观测变量4是共同因子3是独特因子， A 是 
因子负载。因子负载是指共同因子一个单位的变化所引起 
的观测变量的变化。 

通常，观测变量和潜变量之间的关系可以用下面的矩阵 
方程来表示(《验证性因子分析》中的方程 2.4) : 


x = A^+5 [1.2] 

x 是观测变量向量4是共同因子向量，5是独特因子向量。 
从统计上来讲，测量模型的任务就是通过方程 1.2 所定义的 
观测变量和潜变量之间的关系来解释观测变量间的相互关 
系，这个模型是《验证性因子分析》一书的主题。 

结构方程模型已被广泛应用于社会科学和行为科学。 
其最简单的形式就是只含有一个因变量的多元回归模型。 
例如在方程 1.3 中,所有的变量都假定已经均值对中。 

y + P2OC2 +e [1.3] 

^是因变量，: r 是自变量，和 V 的关系是由系数怂和仏来 
表示的^是误差项，表明 ：r 不能完美地预测在第3章和 
第4章中，方程 1.3 将用以下方程来 标示： 
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7l ~7u6l 7l2?2 + Cl 


图 1.2 协方差结构模型中的结构部分 

我们也可以通过包含递归或非递归因果关系的联立方 
程系统来构建更符合现实的模型。例如图 1.2 所展示的模 
型。该模型表明观测变量是由观测变量和^决定 
的，方程中的误差项 G 表明这三个变量并不能完全地解释 
%。类似地，％是由观测变量％、和$3决定的，同时包含 
误差项匕。这些结构关系可以通过下列方程来 表示： 

7 】 = /?12^2 + yilfl + 7l2^2 Cl 

rj 2 —^2\ rj \ + 722$2 723^3 + ^2 [1.4] 

或者更为普遍的结构方程 形式： 

rj = Br 1 + r ^+； [1.5] 

I ?指不含测量误差的观测因变量的向量集 j 是不含测量误 
差的观测自变量的向量集， f 是误差向量集， B 是因变量之间 
的回归系数矩阵， r 是自变量与因变量之间的回归系数矩阵。 
多元回归模型、路径分析、联立方程模型以及追踪数据分析 
等模型都是方程 1.5 的特例。我们将在第3章中详细介绍这 
些模型。 
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因子分析和结构方程模型是相互补充的。这是因为结 
构方程模型关于观测变量没有测量误差的假定通常是不现 
实的，所以需要引入观测变量的测量误差。而运用因子模型 
的研究者通常会对潜在因子的结构关系感兴趣，所以需要将 
结构方程模型和因子分析模型联系起来。协方差结构分析 
就是为了满足这种需要而产生的，它同时考虑了因子模型中 
变量的测量误差和结构方程中的误差。 



图 1.3 协方差结构模型中的测量部分与结构部分 


图 1.3 是协方差结构分析的一个简单例子，图中的结构 
关系和图 1.2 相同。不同的是协方差结构分析中的结构模型 
所表示的是潜变量间的结构关系，而非观测变量间的结构关 
系。在协方差结构模型中，潜变量和观测变量的连接方式和 
因子模型是一样的，即通过带有箭头的线条将方框中的: r 和 
: V 与圆圈中的 7 和 f 连接起来。其数学表达与方程 1.1 是相 
同的，对于因变量而言，其测量模型如下： 


y \ = A 'n 71 + f 1 

^3 = A 32 7^2 + (3 


yz = 又 & 7 1 + f2 

3^4 = A 42 +(A 
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对于自变量而言，其测量模型如下： 

•2^1 —An $1 + 5i oci — Aji$i ~l~ Sz 
工3 = 又二2?2+汐 3 工4 = 又:2令2+汐4 
工5 =又【3冬3 +汐5 工6 =义: 3 冬3 +占6 

用更普遍的形式来表达，协方差结构模型包含三个方 
程。第一，用结构方程来设定潜变量之间的因果 关系： 

T] =B?| + + ^ [1.6] 

这和方程 1.5 的形式一样，不同的是这里的 ry 和 《所 指的是 
潜变量，而非观测变量。 

第二个和第三个方程和方程 1.2 的形式相同，是一对测 
量模型。第一个测量模型通过因子负载矩阵&将观测变量 
x 和潜变量6连接起来： 


x = A,^ +8 [1.7] 

5表示 x 的测量误差。第二个测量模型通过因子负载矩阵 
A, 将观测变量^和潜变量 I ?连接起来： 

y = A x r| + € [1.8] 

€表示 y 的测量误差。测量方程 1.7 和 1.8 连同结构方程 1.6 
共同组成了我们将在第4章中讨论的协方差结构模型。 




: 












m 




c 


i 
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协方差结构模型的测量部分是由一对验证性因子分析 
模型构成的。本章将会介绍两个因子模型的具体形式及相 
关假设。由于参数的设定、估计、检验和解释都与《验证性因 
子分析》一书中的因子模型是一致的，因此，本书将不再讨论 
这些问题。 



第 2 章测量模型 


第 1 节 I 测量模型的设定 


协方差结构模型包含两组因子模型。纟包含 X 中 q 个观 
测变量的 s 个共同 因子; ly 包含 y 中 p 个观测变量的 r 个共 
同因子。在协方差结构模型中， f 指结构模型中的外生潜变 
量(我们将会在第3章中定义相关术语）。观测变量和共同 
因子通过下列因子方程连接 起来： 


Ar^ +5 

[2.1] 

A v rj +€ 

[2.2] 


A , 是: r 在$上的因子负载矩阵 (qX s )， 其中 Afj 表示观测变 
量^在$上的因子负载，5是 I 变量的独特因子（或测量误 
差）的向量集合 ( qXl )。 同样，人，是^在7上的因子负载，其 

中 M 表示观测变量 J ，在％ 上的因子负载—是^变量的独特 
因子(或测量误差）的向量集合 (PX 1 )。 所有这些变量都假 
定已经均值 对中： 

E(x)=E(5)=0 E(^) =0 
E(y) =EU) =0 E(r/) =0 


每个方程中的公因子和独特因子都假定是不相关的，具体 
来讲： 
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E (鉍 '）=0 或者 E (5^)=0 
E (? jf ') =0 或者 E(eq r ) =0 

表 2.1 简要列出了相关假定。协方差矩阵 OXsXs ) 包含 
了 6的方差和协方差，协方差矩阵 找 (qX q ) 包含了 5的方差 
和协方差，协方差矩阵0 (pX p ) 包含了 f 的方差和协方差， 
7的方差和协方差则包含在对称矩阵 COV ( ry)( r X r ) 中。我 
们没有给矩阵 C () V ( T |) —个独特的字母(如0代表 f 的方差 
和协方差矩阵)是因为在协方差结构模型中， 00 X (1^) 是通过 
其他参数来定义的，这点我们会在第3章中详细介绍。 


表 2.1 协方差结构模型测量部分 


矩阵 

维度 

均值 


协方差 

维度 

注 释 


(sX 1) 

0 

0 = 

e(40 

(s X s) 

外生共同因子 

X 

(qX 1) 

0 

2^, 

= E(xx ’） 

(q X q) 

外生观测变量 

Xr 

(q X s) 

— 


— 

—— 

乂在€上的因子负载 

8 

(qX 1) 

0 

©5 = 

= E(S 8 , ) 

(qXq) 

x 的独特因子 

V 

(rX 1) 

0 

COV(ij) = E(rjT]') 

(rX r) 

内生共同子 

y 

(pX 1) 

0 

2I v> : 

= E(yy ’） 

(pX p) 

内生观测变量 

A, 

(p X r) 

— 


— 

—— 

y 在 |? 上的因子负载 

€ 

(p X 1) 

0 

0 f = 

: E(^) 

(p X p) 

y 的独特因子 


注 :因子方程: 


协方差 方程: 


x = 4-6 

y = A . v ij+f 



AvC()V(7j)A ； -h0, 

A v COV(ij, §)Ax 

A., C()V(? , rj)A ； 



[ 2 . 1 ] 

[ 2 . 2 ] 


[ 2 . 3 ] 


假设： 

(1) 所有变量的均值为 0( 对 中）： E ( x ) = E (6) = 0； E (4) = 0； E ( y )= 
E ( f ) = 0； E ( i |) = 0： 

(2) 共同因子与独特因子不 相关： E (鉍 '）=0 或 E (6^) = 0; E (^ ； ) = 
0 或者 E ( erj ')=0 ; E (》'）= 0 或 E (式 ’）= 0; E ( i ^) =0 或者 E (6 V ) =0; 

(3) 不同方程中的独特因子不 相关： E ( ed r ) = 0或 E (€^) = 0 o 
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协方差结构模型测量部分的设定包含在1， A v ， 4>， 
0,，0, 等五个参数矩阵中，可以限定某些参数相等或者限定 
某套参数相等。 

如果将方程 2.1 中的岛替换为 A ， 方程 2.2 中的 A v ， r / 
和 COV ( iy ) 替换为 A ,， A ， 《， O , 5,那么方程 2.1 和方程 2.2 
中的每一个因子模型都和《验证性因子分析》中所描述的模 
型是相同的。但是，将两个因子模型连接起来需要一些新的 
假定。 



协方差结构 模型: IXSREL 导论 


第2节 I 两个因子模型之间的关系 


尽管有些变量可以在两个因子模型之间是相关的，然而 
有些假定则是不能相关的。观测变量 x 和3^可以相关，它们 
的协方差包含在矩阵 UqXp ) 中，其中第 ( i ， j ) 个元素表示 
之间的协方差。类似地，外生潜变量 f 和内生潜变 
量7/也可以相关，它们的协方差则包含在矩阵 COV ( ? , $) 
(rX s) 中 。 

cov(6 2 , 7 3 ) 



图 2.1 协方差结构模型中的测置部分 


我们假定在同一因子模型中，独特因子与共同因子是不 
相关的，在不同因子模型中，独特因子与共同因子也是不相 
关的，即 
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E (私 '）=0 或 E ( g ')=0 
=0 或 E(5i/ ) =0 

尽管5与< 矩阵的内部元素是可以相关的，但是8和 （ 之间是 
不能相关的，即 


E(&，）=0 或 £(€5，）=0 

最后假定观测变量: T 的因子负载不能落在潜变量 7 上，而观 
测变量3^的因子负载也不能落在潜变量$上。 

图 2.1 展示了测量模型的相关假定。圆圈代表潜变量， 
矩形代表观测变量，连接两个变量的双箭头表示这两个变量 
是相关的，单向箭头表示起点变量影响箭头指向的变量。中 
间较粗的实线将两个因子模型分开，观测变量间的相关或者 
因子负载不能跨越这条线。 
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协方差结构 模型: USREL 导论 


第 3 节 I 协方差结构 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 


既然彳和 f 都是不可观测的，模型参数就必须通过观测 
变量的方差和协方差来估计。验证性因子模型的估计方法 
(详参《验证性因子分析》)有助于我们考虑如何将观测变量 
之间的方差与协方差和模型参数联系起来。 

定义5：为包含了总体中观测变量方差与协方差的矩阵： 


Z = E 



"_y " 


丫 

t 


E(yy ) 

E(yx )" 


_ X _ 


_ X _ 



_ E ( x /) 

E ( xx ')_ 


用方程 2.1 和方程 2.2 替换其中的 x 和 y: 



E [( Ayty +€) (+£)’] 

E [( A v 7 y +6)( Aj +6 )l 

— E [( A ,§ + 5)( A v ?y + €)'] 

E [( A ^+5)( A ^+6) / ]_ 


运用模型中关于某些变量不能相关的假定，然后通过一系列 
的数学变换(将此矩阵转置、相乘，然后取其期望），矩阵2：可 
以被写作如下 形式： 



A v COV ( tj ) A v + 0. 

A .', COV (”， ^) A 7 

l A , COV (^ rj)Ay 

A r OA r +0,. 


[2.3] 


读者或许想知道方程 2.3 的推导过程•但是其推导过程 
对于理解这个模型并非必需的。最为重要的是方程 2.3 将观 


第 2 章测量横型 


测变量: r 和 的方差与协方差分解成了因子负载矩阵&和 
A _、 、共同因子 f 和7的方差与协方差矩，以及独特因子8和 c 
的方差与协方差矩阵。协方差矩阵2：会受到这些参数矩阵 
的相关假定的影响。而模型参数的估计就是通过方程 2.3 再 
现样本中的方差与协方差矩阵。 
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协方差结构 模型： IASREL 导论 


第4节 | 小结 


本章介绍了协方差结构模型中的测量部分。尽管验证 
性因子分析模型在《验证性因子分析》中是一种强大的统计 
模型，但其在本章中只是用于将观测变量与潜在因子联系起 
来，其目的是设定这些因子之间的关系。下一章将介绍协方 
差结构分析中的结构部分。 


; 


结构方程模型 
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协方差结构 模型: LISREL 导论 


结构方程模型负责设定协方差结构分析中潜变量的因 
果关系。为简便教学，本章假定协方差结构中的潜变量都是 
可以直接观测的。这种简化有助于读者理解协方差结构分 
析中的结构部分，同时读者也可以借助计量经济学文献中关 
于观测变量结构方程模型的相关知识（按照难度递增的顺 
序，参考 Hanushek &• Jackson , 1977； Wonnacott ， 1979； 

Theil , 1971)。 本章关于结构方程模型的设定、识别及解释 
的讨论，都可以直接应用于第4章协方差结构分析中的结构 
方程模型。 

假定不含测量误差的第二个优势是很多协方差结构分 
析软件在估计含有参数限定的结构方程模型时极为便利。 
结构方程分析中关于参数的限定在社会科学中很常见，尤其 
是在追踪数据模型中。确实，目前为止，许多学者 (Kessler & 
Greenberg , 1981) 建议追踪数据模型可以通过协方差结构模 
型软件来估计，但却没有给出具体的操作方法，本章的第二 
个目的就是填补这一空白。 



第 3 牽结构方程模型 


第1节 | 数学模型 


结构方程模型设定了一组变量间的因果关系。模型中 
被解释的变量称做内生变量，内生变量可以由其他内生变量 
或者外生变量来解释。外生变量则是由模型之外的变量决 
定的，不能被模型解释。 

定义1/为 (rX 1) 维的内生向量集合4为 (sX 1) 维的外 
生向量集合。模型假定这些变量是通过线性结构方程系统 
联系的。 


iy = B 7 y + r^+C [3.1] 

B 是一个 （rX r ) 维内生变量关系的系数 矩阵； r 是 （rX s ) 
维的外生变量与内生变量关系的系数矩阵4是 （rX 1) 维 
的用来表明内生变量不能被结构方程所解释的误差向量 
集合。 

方程 3.1 被视为结构方程，是因为它描述了变量间所假 
定的因果结构。限定 B 和 r 中的部分元素为0则意味着相 
对应的变量间没有因果关系。限定 r 矩阵中的第 ij 个元素 
为0(>^=0)则表明外生 变量& 对内生变量 7 i 没有影响。相 
似地，如果 B 矩阵中的第 ij 个元素为 0 (ft = 0)， 则表明内生 
变量 7 ，不受内生 变量％ 的影响。 B 矩阵的对角线为0,表明 
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协方差结构 樓型: LISREL 导论 


内生变量不受其自身的影响。 

方程 3.1 两边同时减去 Bry ， 得到 i | — Bi 7 =rg + f ， 将 B 
定义为 (1— B )， 则： 


Brj = T ^ + ； [3.2] 

尽管方程 3.2 在计量经济学的文献中更为常见，可是方程 3.1 
更便于解释。 B 的符号为正则表明两个内生变量之间存在 
正向关系，可是6的符号为正则表明了负向关系。由于6更 
便于表示一组结果，因此通常用来代表 ( I 一 B ) 。 

本章中的结构方程模型假定所有变量的均值为0,即 
E ( ry ) = E ( f ) =0, E (4) =0。由于结构参数包含在 B 和 r 矩 
阵中，并不会受到这一假定的影响，所以这个假定并不影响 
结构方程模型的普适性。 

正如在因子模型中假定独特因子和共同因子不相关，结 
构方程模型中同样假定内生变量方程的误差和内生变量自 
身也是不相关的，即 E (私 '） =0或 E ( d ) =0。 

结构方程模型还假定 6 = (1 — B ) 是非奇异矩阵（即6 1 
存在），但是这个假定要求并不严格，仅是表明模型中没有冗 
余的方程。 

表 3.1 列出了结构方程模型的相关矩阵。残差^的协方 
差包含在 (rX r ) 维的对阵矩阵 f 中，既然 f 假定均值为0,即 
E (^) =0,矩阵平可以被定义为 ¥ = E (« f )。 尽管除对角线 
的所有元素都可以限定为0以表示两个方程中的误差是不 
相关的，但矩阵平的值是不知道的。 OS ( sXs ) 维的对称矩 
阵，表示外生变量的协方差。同样由于假定外生变量的均值 
为0, 0> = E (茲'）。 
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表 3.1 协方差结构模型中的结构部分 


矩阵 

维度 

均值 

协方差 

维度 

注 释 

V 

(rX 1) 

0 

COV ( tj ) = ECtjtj) 

( rXr ) 

内生变量 

4 

(sX 1) 

0 

0= E ( 妓 ’） 

(s X s ) 

外生变量 


(rX 1) 

0 


(rX r ) 

方程误差 

B 

(rX r ) 

— 

—— 

—— 

内生变量间的直接效应 

•• 

B 

(rX r ) 

— 

—— 

— 

( I - B ) 

r 

(r X s ) 

_ 

— 

— 

外生变量对内生变量 
的直接效应 


注 :结构方程： 





1 ； = Brj + R + < 

[3.1] 




[3.2] 

简化 形式： 

7 j = B ' r 4 -hB 1 

? 

[3.4] 

协方差 方程： 




2 = 

uror '+ v ) S — 1 2 

b ] ro 一\ 

[3.5] 

L or ' B ’ 1 

rv - 」 

假设： 





(1) 所有的变量都假定均值为0: E ( iy ) = E (^) = E ( f ) = 0； 

(2) 内生变量的误差和内生变量自身也是不相 关的： E (? f ') = 0或 
E(f^) =0 ； 

c . 没有冗余的方程： (1 — B ) - 1 存在。 

通过这些假定，我们可以定义一系列的协方差矩阵。方 
程的误差项？包含在对阵矩阵 f (r X r ) 中，由于假定误差项 
的均值为0,即 E ($)=0， W 可以写为 f = E ( W )， 尽管假定 
平矩阵中的非对角线元素为 0( 即不同方程中的误差项不相 
关）， W 的值仍是未知的。外生变量的协方差矩阵包含在对 
阵矩阵 O ( sXs ) 中。 

下面我们通过几个例子来展现这些假定及结构方程模 


型的灵活性。 
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协方差结构 模型: LISREL 导论 


例 3.1: 多元回归 


下面的方程表示含有一个内生变量和三个外生变量的 
结构方程 模型： 




[7i] = [Q] [7 丨 ]+ H/ny^yi.sH 62 H - [?i] 

_^3_ 


或者写为： rj \ =7 h 6 i + yn ^2 + yi 3 f:i + ? i ( 这里及本书的其他 
地方，含有下划线的数字表示其对应的参数等于该数字 ）。 o 
矩阵包含了总体中 f ,， f 2 的方差与协 方差： 


■VAR(^) COV(f M e 2 ) COVC6,, $ 3 )" 

0= COV(6 2 , ft) VAR ( 匕） COV(f 2 , e 3 ) 

COV(e 3 , 6i) COV($ 3 ， f 2 ) VAR ( 匕 ） _ 

+ 11 《12 多 13 

= j>2\ (j>2Z 彡 23 

jj>Z\ ^32 於 33 _ 

平仅包含方程中误差项的方差:平 = [%]。例 3.1 是最简单的 
结构方程模型，等同于只有一个方程的多元回归模型。// [1] 


例 3.2: 追踪数据模型 


通过增加一系列的方程，例 3.1 中的模型可以变得更加复 
杂。例 3.2 来自惠顿 ( Wheaton ，1978) 关于心理失调的追踪数 
据分析模型的变形。图 3.1 设定了父亲的社会经济地位 (&) ， 
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研究对象三个时期的社会经济地位（％， 72 和％)，以及研究 
对象在两个时点上的心理失调症状(％和％)。原文中有两个 
关于心理失调的指标，本例中只有一个指标，即假定该变量没 
有测量误差。第4章会介绍含有两个测量指标的原模型。 



注: FSES 为父亲的社会经济 地位; SES ] 为研究对象吋点1时的社会经 
济 地位; SES 2 为研究对象时点2时的社会经济 地位; SES 3 为研究对象时点 
3时的社会经济 地位; PSY 2 为时点2时心理失调症状的 数目； PSY 3 为时点 
3时心理失调症状的数 H ; PSY 2 为时点2时精神生理失凋症状的 数目； 
PSY 3 为时点3时精神生理失凋症状的数目。 

资料 来源: Wheaton , 1978。 

图 3.1 心理失调的结构模型 


图 3.1 中的结构方程模 型为: 


7'~ 


"0 0 0 0 0 ' 


V 


>n' 





/? 21 0 0 0 0 


72 


y?i 



V' 


^31 /?32 0 0 0 



+ 

y 3 i 


63 

V 、 


/?41 ^42 /^43 0 P 


74 


y.n 



-7 - 


fis\ /?52 /^53 0 0 _ 


-75- 


_/ 51 _ 


^5 - 


[3.3] 
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协方差结构 模型: LKREL 导论 


既然只有一个外生变量 ,0 = [ VAR (^)]。 

方程中残差的协方差包含在矩阵¥中。此处平可以有 
两种设定 :第一 ，假定方程中所有变量的误差都不相关，即将 
f 限定为对角线矩阵。第二，假定其中某些残差是相关的， 
这更符合实际情况，在图 3.1 中用虚线的双箭头来表示。这 
些设定表明在同一时期内，研究对象的社会经济地位与其心 
理失调的症状是相关的。如果模型没有正确设定同时影响 
社会经济地位和心理失调症状的变量，就会发生这种情况。 
由于此种错误模型设定所造成的误差也会出现在时点2的 
G 和 G 中以及和时点3的匕和&，因此这两对误差也是相 
关的。我们会在模型的设定和估计中讨论¥矩阵设定的重 
要性。含有这些限制少矩阵的设定 如下: // 

0 „ 0 0 0 0 " 

0 ip 22 023 0 0 

平= 0 032 033 0 0 

0 0 0 044 045 

_ 5 5 9 . 054 055 _ 

例 3.3: 同时因果关系 

本例来自经常被引用的邓肯及其同事 （Duncan et al .， 
1971) 的模型。在这个模型中（图3.2)，研究对象及其朋友的 
特征被用来预测本人及朋友的期望。原文中用了两个指标 
来测量期 望:职 业期望和教育期望。本章只采用教育期望这 
一个指标，％为研究对象的教育期望，％为研究对象的朋友 
的教育期望。模型中共有六个外生变 量:匕 测量研究对象的 
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父母对其的期 望，^ 测量研究对象的智力水平 （IQ),f 3 测量 
研究对象的社会经济地位 (SES)，$ 6 测量研究对象朋友的父 
母的期望，$ 5 测量研究对象朋友的智力水平 （IQ)， 匕 测量研 
究对象朋友的社会经济地位 (SES)。 



注: R-PASP 为研究对象父母的 期望; R-IQ 为研究对象的 IQ ; R-SES 为 
研究对象的社会经济地位; F-SES 为研究对象朋友的社会经济 地位； F-IQ 为 
研究对象朋友的 IQ ; F-PASP 为研究对象朋友父母的 期望; R-ASP 为研究对 
象的期望; F-ASP 为研究对象朋友的期望。 

资料来源 : Duncan et aL , 1971。 

图 3.2 同辈对期望影响的结构方程模型 


这个模型的结构方程表达 如下: 
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即研究对象的期望(％)会受到其父母的期望 (匕） 、智力水平 
( 6 2 )、社会经济地位（匕）及其朋友的社会经济地位 （匕） 的影 
响 ( ZiifO ， y 12 /0， / 13 ^ 0, y 14 #0， y l5 = 0, yi 6 = 0)。 同 
样，研究对象朋友的期望 （ 72 ) 会受到其父母的期望（匕）、 
智力水平（匕）、社会经济地位（匕）及其研究对象的社会经 
济地位的影响 （ f 3 )(y 23 关 0, / 24 ^ 0, y 2S ^ 0, y 26 关 0, 
y 21 = 0, y 22 = 0)。该模型还设定研究对象及其朋友的期望 
相互影响(/? ]2 关0, /? 2] #0),并预测研究对象期望的方程的 
误差与预测研究对象朋友期望的方程的误差是相关的，这是 
因为模型可能没有包含同时影响研究对象及其朋友期望的 
变量。// 
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第2节 | 协方差结构 


我们可以通过方程 3.1 来定义内生变量之间的协方差、 
外生变量的协方差 （中） ，及残差项的协方差（少)。由于假定 
B 是非奇异矩阵，因此方程 3.2 两边同时乘以6的逆矩阵，得 
到方程 3.4: 


t]=B + B [3.4] 

方程 3.4 用外生变量和方程误差来解释内生变量，是结构方 
程模型的简化形式。 

既然假定 E(r,) =0,那么1；的协方差矩阵就等于 
ECttV )， 通过简化方程3.4, 7?的协方差可表 达为： 

COV ( iy ) =E(iyry ， )=E[(fi- 1 r^ + B- , p(B- , r^+B 
= E[(6— 1 B ) + (B" 1 B ) + 

(ft— 1 g’r’6’ _1 ) + (B 

去掉上面的数学期望，并通过 £(^)=0 及关于4> 和少矩阵 
的定义，得到如下 形式： 

COV ( rj ) = B ror fi -1 +B 

=b 1 (ror+^)B , ~ 1 

因此内生变量的协方差可以通过结构参数 r 和6、外生变量 
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协方差结构 模型: USRKL 导论 


的协方差，及误差的协方差来表达。 

可以通过相似的方式得到6和 7 的协方差矩阵。因为 
E(f)=0, E(iy)=0, 所以 COV(i;，g) =E(r^') ，用简化形式 
B r^ + B ^ 代替 q ，通过 4>= e ( g ') =0 的假定 得到： 

cov ( i |, = E (?；^) = e [( b - 1 i ^ + 

= E[fi-irW + S-W] 

=b 1 ro 


通过上述结果，我们可以定义包含 f 和 7 的方差和协方 
差的矩阵 2 Xr + sX r + s ): 



"COV(iy) 

COV(f|, §)" 


COV(^, r\) 

COV(^) 


6 (ror +^)S 一 1 

frro 

or’ft’ 

-1 

o 


[3.5] 


2：矩阵是通过结构参数矩阵 B 和 r 、 协方差矩阵 cp 和¥来定 
义的。结构方程模型可以限定结构参数以及协方差矩阵 z 
中某些元素的值。在实际中，我们并不知道总体中2：的值， 
但是样本的协方差矩阵 S 是已知的。结构方程估计的过程 
就是通过方程 3.5 找到一个与样本的协方差最接近的矩阵。 
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第3节 | 结构方程模型的类型 


结构方程模型的识别和估计主要取决于 B 和平 矩阵的 
形式。接下来，我们首先区分三种形式的 B 矩阵和两种形式 
的少矩阵。 

B 矩阵的三种类型 ：（1) B 为对角线矩阵，内生变量只受 
外生变量的影响，内生变量之间没有关系； （2) B 为三角 
矩阵，表示内生变量之间存在因果关系，但是这种关系是单 
向的，如果％影响％，则％不会影响，即如果 ft # 0,则 
ft =0;(3) B 矩阵在对角线外有不受限制的元素，这说明内生 
变量之间可以相互影响，也就是说，％可以影响 7i ， ％也可 
以影响％，即如果则 ft 可以等于0,也可以不等于0。 

平矩 阵的两种类型：（1)平为对角线矩阵，这表明所有的 
误差都是不相关的，即只要 i 不等于】，则少, i =0; (2)^ 
为对称的非对角线矩阵，即至少有两个方程的误差项是相关 
的。也即，至少有一个少6=%,且不必然等于0。 

根据 B 和平矩阵的组合，共有六类结构方程模型。 

第一类模型是最简单的情况，即 B 和平都是对角线矩阵。 
如果只有一个方程，就是多元回归模型(如例 3.1 )。如果多余 
一个方程，则方程之间是相互独立的，即某一方程的估计与其 
他方程无关，每个方程都可被看做单个的多元回归模型。 
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第二类模型中内生变量之间不相关•但是方程的误差项 
相关 （ B 是对角线矩阵，少矩阵中至少有一个对角线外的元 
素不为0)。如果两个方程的外生变量完全不相同•就是这种 
类型•即“看似不相关的回归”(参见 Kmenta . 1971)。 

第三类模型中 B 为三角矩阵•表明内生变量可以影响其 
他内生变量，但是内生变量之间不存在相互因果关系，即如 
果％影响则％不会影响 7 i ( Ai =0) 。 当这种形式 

的 B 矩阵与对角线平矩阵组合在一起时，就组成了递归模 

1 % 

型。例如 3.2 中的 B 为三角矩阵，当设定少为对角线矩阵 
时，这个模型是递归模型。 

第四类模型与递归模型不同的是， V 为非对角线矩阵。 
这种情况通常发生在追踪数据模型中，即不同变量的误差在 
同一期内相关，或者是同一方程的误差项在不同时期内相 
关。在例 3.2 中， B 为三角矩阵，如果设定¥为非对角线矩 
阵，这个模型就不再是递归模型。 

第五类模型是，尽管在同时因果模型中 （ B 矩阵不是三 
角矩阵或者对角线矩阵)通常假定不同方程的误差项是相关 
的（即 W 为非对角线矩阵），但这并非必需的。如果能够证 
明，无论 B 矩阵是什么形式，方程中的误差项都是不相关的 
(即少为对角线矩阵），就是这种情况。 

第六类模型是，当内生变量之间存在同时因果关系并且 
不同方程的误差项是相关的。这类模型在计量经济学中被 
称做非递归同时因果模型。例 3.3 就是这类模型。 

请记住这些类型的模型，在模型的设定和估计时可能用 
到这里的讨论。 
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第4节 | 模型的识别 




wmmmmm 


如果结构方程模型不能识别，则意味着模型参数有多个 
解可以再现观测数据。在这种情况下，研究者无法从中选 
择，因为所有解都是有效的或者都是无效的。模型的识别可 
以通过限定 b , r ， 平等矩阵的参数来实现。所谓模型是可 
识別的，是指在既定参数设定下，有且只有一组解符合模型 
的要求。 

无限定的结构方程模型（比如对 b ， r ， 平矩阵没有任何 
限制)是不能被识别的（证明请参考 Wonnacott Wonnacott , 
1979： 462)。当限定 B 和屮矩阵为特定形式时，模型通常是 
可以识别的。具体来说，不相关的多元回归（第一种类型）、 
看似不相关的多元回归（第二种类型）以及递归模型（第三种 
类型)都是可以识別的。除非对 b ， r 和平矩阵有足够的限 
定，否则后面三种类型的结构方程模型是不可识别的。接下 
来，我们介绍什么样的限定是可行的，以及如何证明模型是 
可以识别的。 

联立方程模型(第六种类型）的识别可以通过设定 B 和 r 
矩阵中的部分元素为0来实现。如果限制 B 矩阵的某一元 
素为0,比如 ft =0, 其意思是内生变量不会影响内生变量 
rj ' 0 如果限制 r 矩阵中某一元素为 0, 比如 yy^O, 其意思是 
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外生变量$不会影响外生变量6。费希尔 （ Fisher ) 在1966 
年就已经讨论了此类模型识别的限定条件。 

最简单的情况是测量阶条件 （order condition ) 0 阶条件 
是指如果系统中的某个方程是可以识别的，那么该方程中已 
知参数的数目一定大于或等于方程组的数目减去1。换言 
之.在 b 和 r 矩阵中，在给定行中限定为0的系数一定大于 
或等于方程组的数目减去1。这是一个必要条件，如果违反 
了该条件，方程一定是不能被识别的；而即使满足了这一条 
件，方程依旧有可能无法识別。下面我们通过例子来阐明这 
一条件。 

例 3.3: 阶条件 

邓肯 ( Ducan ) 、哈勒 （ Haller ) 和波特斯 （ Portes ) 设立了一 
个典型的联立方程组模型。结构系数的矩 阵为： 

_/^21 P _ 

和 

7n 712 713 7l4 0 0 " 

r — 

_ 0 0 723 y 725 726 _ 

r 矩阵的第一行设定 y 15 = y 16 = o , 这说明模型假定外生变量 
^和^对内生变量 71 没有影响。因此％的方程中已知参 
数的数目大于方程组的数目减去 1(2 > 1)。所以该模型满 
足了阶条件。同样，&的方程也满足了该条件。// 

秩条件 (rank condition ) 是结构方程模型可识别的充分 
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必要条件。定义为排除6中该方程所包含的行•然后删 
除在该行中所有不为0的列所形成的矩阵。再以相同的方 
法获得矩阵 r 。 秩条件就是当且仅当矩阵 [ ij * | r s ] 的秩等 
于方程组的数目减1时 [3] •该方程组才是可识别的。虽然在 
实际中我们经常忽略这个条件，但这种忽略是很不明智的。 
尽管由于6和 r 是未知的(它们是在模型可识别的情况下估 
计出来的），但有时我们是可以知道当[以 | 广]的秩小于方 
程组的数目减1时，给定的方程是不可识别的。 

例 3.3: 秩条件 

为了检验内生变量％的方程是否满足秩条件，首先通 
过删除6和 r 矩阵中的相关元素，构建矩阵和 r % 

•• ^ 

= 

「I 1 _ 


和 


b b 

r = - 

N N y 25 y 26 _ 

相应地， [6^ r ] 等于 [ y 25 , 726] o 当总体中的 y 25 和 y 26 同时 
等于 o 时， [6# |广]的秩等于1。尽管不是很确定，我们还是 
可以推 导出％ 的方程是可识别的。同样，秩条件也可以应 
用到 72 的方程中。// 

例 3.2: 秩条件和阶条件 


惠顿 ( Wheaton ，1978) 的模型（图 3.1 和方程 3.3) 是个典 
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型的考虑到同一时期内不同方程的误差项存在相关的追踪 
数据模型。如前文所述 . B 为三角矩阵，如果假定平为对角 
线矩阵，则该模型为递归 模型. 是可识 别的； 如果假定少为对 
称的非对角线矩阵，尽管 B 为三角矩阵，该模型也不是递归 
模型，而且此模型也未必是可以识别的。我们可以应用阶条 
件来检查该模型是否可以识別。 


因变量 已知参数 


已知参数是否> 
方程数目减1 


方程数目减1 


7 i 4 是 

^2 3 否 

7 s 2 否 

7, 1 否 

75 1 否 



由上可见， 从&到 W 的方程都是不可识别的。可以应用秩 
条件来检验的方程是否可以识别。首先构建 [w ir ] 
矩阵： 


[护 I r]= 




§ 





1 

0 

0 

0 



卢32 

1 

0 

0 



P \2 

/?43 

1 

0 



/?52 

卢53 

0 

1 



1 

々32 

^52 


0 


1 




0 0 
0 0 
1 0 
0 1 


[ B « | P ] 的秩为4,因此，当不限定矩阵¥时，的方程是可 
以识别的。然而如果限定少矩阵，剩余的方程或许是可以识 
别的，但是我们不能通过阶条件和秩条件来检验。 

除了限定参数为0夕卜，等同性限定(如 y 23 = y 25 ) 、非线性 
限定及方程误差项协方差限定（如少 12 =少 21 =0) 等，也能够 
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使模型变得可以识别。当运用上述限定时，有必要通过观测 
变量的方差和协方差来证明模型参数是可解的。下面我们 
通过例子来展示如何去证明模型是可以识别的。 

例 3.2: 模型识别 

在惠顿 （ Wheaton ，1978) 的模型中（图3.1)，如果不限定 
f 矩阵是很有问题的。在原模型中•惠顿假定同一时期内的 
方程的误差项是相关的。中矩阵为： 

0 ,, 0 0 0 0 ' 

^ 022 023 0 0 

平= 0 032 033 P 9 

0 0 0 ^14 045 

_ 0 0 0 054 055 _ 

通过限定 W 矩阵中的元素，该模型中的其他方程或许也是可 
以识别的。为了检验模型是否可以识别， b ， r 和平矩阵中 
的参数必须可以通过外生变量和内生变量的协方差得出。 
结构方程组可以写做如下 形式： 


= 

= 

7 \\^\ + Ci 

[3.6] 

72 = 

=馬1 + 

721 6 l + 

[3.7] 

7-3 = 夕31 + 032 7 2 + 

/3 lfl + ^3 

[3.8] 

V ' = 

= 夕41 % + /^2 % +仏3 

73 + 7 u 6 l + 

[3.9] 

75 = 

+ 053 7 3 + 7 51 & + 

[3.10] 


通过阶条件和秩条件，我们可以知道方程3.6(即 y n ) 是可以 
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识别的。方程 3.6 乘以其自身然后取期望 得到： 

E(”0 VAR(i| 1 )= y ii 2 VAR(^)+2 y 11 COV(； m ^)+11/,, 
=y n 2 VAR(4，) + i^ n 

既然 y u 是可以识别的， VARW,) 是已知的， COV(&， 匕）假 
定为0,因此如是可识别的。 

通过相同的方式，首先方程 3.7 两边同时 乘以& ，然后取 
期望： 

K( > =COV(Ty 2 , <5i) =^ ： iCOV(i|,, ) 

+ r 21 cov (^, ^) + cov (^, ^,) 

=p ll COV(rj ] , t)+ y 2 iCC)V(U) 

既然假定 COVWdd^O , 接下来方程 3.7 乘以％然后取 
期望： 

E(if 2 f|i) =COV(f| 2 » i|i) 

= p 2 l COV ( iy M ffd + hCOV %， iid + COV %, 咿) 

COVC^e .) =0,可以看作乘以方程 3.6 再乘以匕，然后取期 
望: E( 7l ^)=y„E(f 1 ^ 2 ),既然假定 E (乞 f 2 ) 和 

E(6( 2 )=4 12 =0 •因此 E( t 7，^)=0 o 现在两个方程中有两 
个未知 变量： 

COV (,7/2 9 ) ~p2\ COV ( Tji » f 1 ) ~t~ /21 COV(f 1 f f I ) 

COV( Tjz ♦ 71 ) = /?2i C()V( 7^! , ) + V2\ C()V(f, , rj\ ) 

通过上面的方程组可以解出 /? 21 和 y 21 ，因此该方程组是可识 
别的。 

仏 2 的证明可以通过方程 3.7 乘以其自身然后取期望 
得到： 
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COV( 7 2, 72)=^ 21 2 C(3V( ?1 , T^ + y/COV% ， ^)+0 22 

+ 2^21 / 2 iCOV (rji * 61 ) + 2 卢 21 COV( % ， （ 2 ) 
+ 2 y 2 iCOV ( f " ( 2 ) 


既然假定与^有关的协方差都为0,而且方程中除了如 2 都 
是可知或可识别的，因此仏 2 也是可以识别的。 

通过既有的限定，还剩下方程 3.8 是不可识別的。将方 
程 3.8 分别乘 以匕， 71 和 7 2,然后取 期望： 

COVCr/” f 】） =p 31 C()V( 7 l ， &)+ 知 COV( ? 2 ，乞） 

+ y 31 COV(f" e^+COVCCa, fi) 
COV( 々 3 , 7j\ ) =^3iCOV(l^i » 7l ) +/?32 COV( 7^2 » 7l ) 

+ 7 31 COV(f 1 , v 〗） +COV(G，％) 
COV (^ 3 , r/ 2 ) =^3iCOV(^] , r/2 、 + P 泣 CXJVirjz ， r/2) 

+ y 31 COV($ 】， 72 )+C()V(f 3 , 仍） [3.11] 

尽管 COV (匕， ％) 和 COV%，&) 等于0,可是 COV%， ?2 )= 
如并 不等于0。因此，在三个方程组中有四个未知元素(/? 31 、 

戸32、731、少23)，该方程组是不可识别的。 

由于每增加一个方程就会增加一个新的未知参数，因此 
试图寻找第四个方程组来解决问题是行不通的。例如： 

COV( 73 , 73 )= Ai'COVCtjim p+i^COVO^ ，仍） 

+ y 3 i 2 COV(f ， f 〗 ）+ 知 2 + 2 /? 31 y 31 COV( 7 " 
61 ) + 2032 y3i COV( rj2 ♦ 61 ) + 2/?32^3i C( )V( rj \» 
^ 2 ) + 2y 3 i COV(f 1 ^ G) + 2/? 31 COV(7^ ， ) 
+ 2々 32 COV(7 2 ，（ 3 ) 

当增加一个新的方程时，就会增加一个新的未知参数如仏3。 
问题由三个方程解四个未知参数变成了四个方程解五个未 
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知参数。因此，在现有的限定下，尽管方程 3.6 和方程 3.7 是 
可识别的，该模型仍是不可识别的。 

惠顿 ( Wheaton ，1978) 在其模型分析时（尽管他的模型 
包含了测量模型)也遇到了相同的问题。为了使模型可以识 
别，惠顿限定戽 2 ，/?,,和 卢 51 等于0,也就是说内生变量中的三 
条因果路径被假定是不存在的，通过这些限定，可以很容易 
地证明该模型是可识别的。 

现在,方程组 3.11 包含三个方程，可以用来解决三个未 
知参数。以同样的方式可以证明模型中的其他参数也是可 
识别的。// 

关于如何通过限定 f 矩阵来证明模型是可识別的例子， 
可以参阅哈努谢克和杰克逊的论述 （Hminshek &• Jackson , 
1977)。 
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第5节 | 模型的估计 


一旦确定结构方程模型是可识别的，我们就可以开始估 
计程序了。如果模型是不可识别的，那么估计出来的参数也 
是没有意义的。 

和模型的识别一样，模型的估计也依赖于 B 和少矩阵 
的形式。对于不相关的多元回归模型（第一种类型），一般 
最小二乘法是最适合的估计方法。如果 B 为对角线矩阵， 
但少不是对角线矩阵（第二种类型），那么其估计方法取决 
于外生变量。如果在所有的方程组中外生变量都是相同 
的，一般最小二乘法就可以提供无偏、有效的 估计； 如果不 
相同，最小二乘法就可以提供无偏、一致但不是有效的 估计. 
在这种情况下，一般最小二乘法可以提供有效的估汁。 [4] 对 
于递归模型(第三种类型），一般最小二乘法可以提供一致和 
有效的估计，如果模型不包括滞后的内生变量，那么 一般最 
小二乘法的估计还是无偏的。对于没有限定 平矩阵 的联立 
方程组而言，可以采用以下几种估计方 法：二 阶段最小二乘 
法、工具变量、完全信息最大似然估计法及三阶段的最小二 
乘法等。 

计量经济学文献广泛而深入地讨论了每种模型及其 
估计方法，比较简单的介绍参考哈努谢克和杰克逊的著作 
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(Hanushek Jackson ，1977); 更深入的介绍参见克曼塔和 
马林沃德的论著 （ Kmenta ， 1971 ； Malinvaud ， 1970)。这些 
模型可以通过比较常见的统计软件进行估计。对于这些方 
法的介绍不是本书的内容，本书关注那些不能通过常用的统 
计软件进行估计，以及计量经济学文献讨论较少的方法。例 
如，限定¥矩阵中某些非对角线的元素为0或者是相等的模 
型，通常不能通过处理回归模型的软件进行估计，在这种情 
况下，针对协方差结构模型的统计软件提供了便捷的估计 
方法。 

在讨论具体的估计方法之前•有必要区分完全信息和有 
限信息的估计方法(有时候也被称为单一方程方法和系统方 
程方法，参阅 Hanushek & Jackson , 1977)。有限信息估计 

技术在不考虑其他方程的限定下分别估计每个方程。完全 
信息估计技术会考虑到整个系统，并且同时估计系统中的每 
个方程，因此在参数估计时，能够充分利用整个系统提供的 
信息，这在统计上更为有效。另一方面，完全信息估计技术 
也有其局限。既然每个参数的估计都依赖于模型中的其他 
参数•那么它们就更容易受到模型设定的影响。而有限信息 
估计技术尽管不是那么有效，但是由于其分别估计每个方 
程，因此较少受到模型设定的影响。因此，当模型的设定不 
是很明确时，研究者比较偏向有限信息估计技术。 

协方差结构估计软件的程序如下。 [5] 研究者从样本的协 
方差矩阵 S 开始，对角线上的元素为观测变量的方差，对角 
线外的元素为观测变量的协方差。如果数据经过标准化处 
理， S 矩阵就是相关系数矩阵。将 S 看做如下的分割矩阵很 
有 帮助： 
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内生变量之间 
协方差 

内生变量与外生 
变量的协方差 

外生表明与内生 
. 变量的协方差 

外生变量之间的 
协方差 


然后通过 6. r ， 少，少矩阵的估计值（方程 3.5) 来估计总体 
中的协方差矩阵 n : 


/\ 

L = 


COV(T7> 

COVCt;, 

.COV((g, ?|) 

COV(^) 


B _1 (ror 

b 1 ro 

or b 1 

o 


带有帽子符号的参数表明该矩阵包含对总体参数的估计，这 
些参数的估计必须满足相关的模限定。估计的过程就是通 
过找&，士， o 的值来产生一个与样本协方差矩阵 s 最为 
接近的总体协方差的估计矩阵 i 。 其估计过程考虑到了所 
有包含了 6, r ， 平， o 的可能的矩阵。许多矩阵被排除了是 
因为它们没有包含对参数的限定条件。定义厂，平'， 
表示所有考虑到模型限定条件的矩阵，其矩阵集合如下： 


v 铸 _ 

' cov (”） • 

COV (” ， ㊁ ） 

• 一 


2 j 

_ COV (名， ry). 

b -' cr*o r * 

COV (^)" 

, + W* )B -, ~ 1 

] 


如果 s ' 

o*r*B 7-1 

’与 s 非常接近，那么可以认为 

o* 」 

b* , r * , 


是对总体参数比较合理的估计。估计的难点在于如何测量 
两个矩阵 ir 与 s 的“接近”程度，以及如何找到6、 r ， 
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，0> 的值来产生一个与 S 最为接近的 2 T 。 

测量给定的 2 T 与样本矩阵 S 接近程度的函数被称做拟 
合函数， F ( S ; 2 T )。拟合函数表明 5 T 是通过 B •，厂 ，中、 
f * 来界定的，也可以写为 F ( S ; IT ，厂 ， W )。这个拟 
合函数包含了所有满足 B ， r ， 中， f 约束条件的 IT ，厂 . 

^。对于给定的样本协方差矩阵 S ， 那些使得拟合函数 
最小的 B ' 厂，平‘的值就是样本对总体参数的估计， 
记做 S ， f , 分，各。 

协方差结构模型的估计软件中通常运用三种类型的拟 
合函数•它们分別对应于未加权的最小二乘法 （ ULS )、 一般 
最小二乘法 ( GLS ) 和最大似然估计 ( ML )。 

通过未加权的最小二乘法估计方法来估计 B , r ， 4>， 少 
的拟合函 数为： 

F ul < s ( S ； L * )= tr ([( S - E * ) 2 ] 

tr 是一个表明矩阵中对角线元素之和的算法。 ULS 估计方 
法可以在未对观测变量的分布做任何假定的情况下得到一 
致的估计 (Bender 8^ Weeks, 1979 )。 这意味着对于大样本 
而言， ULS 估计是接近无偏估计的。不对观测变量的分布作 
出假定是 ULS 估计方法的优势，但是这种方法也有两个局 
限。第一•没有与 ULS 相关的统计 检验; 第二， ULS 估计的 
值依赖于变量的测量尺度（参考《验证性因子分析》)。 

一般最小二乘法估计方法的拟合函数比较复杂，因为它 
对根据 s - 1 对 S ^ r 之间的差异做了加权处理（详参 
Joreskog &- Goldberger , 1972 )。 GLS 的拟合函 数为： 


F gls ( S ； ^)= tr ([( S -2：- ) S ^ j ] 
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最大似然估计法的拟合函 数为： 

F ml (S ； E.)=tr(Z- _1 S) + [log I 2T I- log | S |]-(r+s) 

bg |2 T I 是 IT 的决定因子的对数。如果 $ 和 T | 为多元正态 
分布， GLS 和 ML 都具有渐进性的特征。 ML 的估计值是最 
接近无偏的.其抽样分布的方差也是最小的，且接近正态分 
布。这意味着如果€和1?为多元正态分布，当样本规模接近 
无限大的情况下， ML 的估计值具有以下特征：（1)样本估计 
值的期望越接近总体参数； （2) 抽样分布的方差会 越小； （3) 
估计值的抽样分布接近正态分布。在协方差结构模型中， 
GLS 估计方法渐进地等同于 ML 估计方法 （ Lee ， 1977； 
Browne , 1974)。两种估计方法都与变量的测量尺度无关， 
且具有统计检验的特性。 

需要注意的是这些渐进性的特征。严格来说，只有当样 
本规模无限大时，它们的特征才能得以保证。实践中常遇到 
的一个问题 是:样 本规模多大才能利用这些渐进性的特征? 
尽管布斯马 （ Boomsma , 1982) 在验证性因子分析中已经取得 
了一些结果(详参《验证性因子分析》)，但对于这个问题，并没 
有确切的答案， GLS 和 ML 估计法都需要正态分布的假定， 
尽管 GLS 对这一假定的要求不如 ML 那么严格 （ Browne ， 
1974)。不幸的是，我们还不是很了解违反这一假定的后果。 


协方差结构 横型: USREL 导论 


第6节 I 实践中的考虑 



没有任何一种估计方法一定优于其他方法。尽管我们 
不需要关注它们是如何在所有可能的 b ， r ， 少，平值中寻找 
使得拟合函数最小的值，但是我们需要注意它们在实际应用 
中的问题。 

第一，寻找局部最小的方法是可行的，这个值使得拟合函 
数比其他最小值更小，尽管这种情况比较少见 (jGreskog & 
Sorbom , 1981 ) 0 

第二，使得拟合函数最小的参数值可能存在于可能值的 
范围之外。例如，估计的方差可能是负值或者相关系数大于 
1，这种情况通常是由于样本规模太小或者不合理的模型设 
定造成的。 

第三，搜索过程可能极为耗时，模型估计的时间取决于 
以下因素：（1)观测变量协方差矩阵中独立元素的 数量; 
(2) 所需估计参数的 数目； （3) 初始值与实际值的接近程度。 
协方差结构估计软件需要对每个要估计的参数设定一个初 
始值。初始值是由研究者提供的猜测值，它们被用来估计第 
一个 IT ，之后软件会重新定义初始值。因此，初始值与实际 
值越接近，估计就会越容易。尽管最新的 LISREL 软件整合 
了产生初始值的算法，但是选择初始值仍然是很困难的。 
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由于 ULS 或者 GLS 的估计值与常用的结构方程模型估 
计方法并不是很一致，因此本章不再做过多的讨论。当满足 
正态分布假定时， ML 的估计与结构方程模型的完全信息最 
大似然估计 ( FIML ) 是一致的(参见 Theil ， 1971)。当违反正 
态分布假定时，完全信息最大似然估计仍然是有效的（参见 
Joreskog &- Sorbom , 1981)。然而，当违反正态分布假定时， 

需要特别留意统计检验。 

最新版本的 LISREL 可以估计初始值。对结构方程模型 
而言，初始值估计类似于二阶段最小二乘法(参见 Joreskog 
Sorbom , 1981)，它在有限信息的情况下，提供了结构方程模 
型参数的一致性估计。 


协方 II 结构 模型： LISREL 导论 


第7节 I 评 估模型拟合度 


本节简要回顾协方差结构模型软件提供的评估模型拟 
合度的相关技术。更详细的讨论可以参考《验证性因子分 
析》。如果结构方程模型是通过传统的软件进行估计的（如 
没有限定等参数相等或者当平为非对角线矩阵时没有施加 
限定），可以运用其他评估模型拟合的方法。关于此类技术 
的讨论可以参见计量经济学的文献。 

检验参数值 

第一步是考虑单个的单数估计值。不合理的估计值意 
味着模型存在问题。通常需要注意以下几点 ：（ i ) 错误的模 
型设定可以导致有偏的参数估计； （2) 由于模型不可识别，而 
产生多个估 计值； （3) 由于样本量太小而导致的估计不精确 
或者渐进性特性无法适用； （4) 输入了不当的相关系数或协 
方差矩阵，或者程序设置错误等。 


估计值的方差 


GLS 和 ML 提供了估计的协方差矩阵。通过这个估计 
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的方差，可以对单个参数进行 Z 检验，以确定参数是否等于 
某个值。这种检验与多元回归模型中的检验方法是一致的。 
需要注意的是，这种检验方法是建立在正态分布的假定上 
的，如果违反了该假定，就需要特别谨慎对待。 

拟合度的卡方检验 

ML 和 GLS 估计法的卡方检验的原假设 （ H 。） 是观测变 
量的协方差矩阵，是由所设定的模型产生的，其备择假设 
( H ,) 是观测变量的协方差矩阵，是一个不受限定的矩阵。拒 
绝原假设意味着模型并没有很好地再现样本的协方差矩阵 
S 0 卡方检验的自由度 为：自 由度 （ df ) =原假设 （ H t) ) 下自由 
参数的个 数一备 择假设 （ H ,) 自由参数的个数。或者更具体 
来说，自由度等于5：中自由参数的数量减去 b ， r ， 少，平中 
自由参数的数量。如果一个模型刚好可以识别，则2：中自由 
参数的数量等于 B ， r ， 4>，少中自由参数的数量。在这种情 
况下，卡方等于0,自由度也等于0,因此卡方检验不能用来 
评估刚好识别的模型的拟合度。 [6] 

卡方差异检验 

正如 F 检验经常用于多元回归中同时检验几个假设 
(Wonnacott Wonnacott , 1979 ： 184—186)，卡方差异的检 

验也可以用于结构方程模型。假设模型 M 1 的卡方值为$， 
自由度为 df , ，模型 M 2 的卡方值为 X 〗，自由度为 df 2 ，且模型 
M , 嵌套于模型 M 2 中，卡方差异 X 2 =$ — X 〗， df ^ df ,— df 2 . 
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协方差结构 模型: L1SREL 导论 


这可以用来检验包含在模型 M 2 中但没包含在模型中的 
参数限定的假设。 

确定系数 

在经典的回归模型中，确定系数是指模型可以解释的因 
变量差异的百分比。在方程中预测 y 的确定系数可以定 
义为： 

VAR ⑹ 

VAR ( 7i ) 

〔，是方程中的误差，在结构方程模型中运用确定系数时需 
要注意，模型所解释的^的差异有可能包含在其他方程中。 

如在例 3.2 中，较大的 R 2 并不一定表明模型的拟合度就 
好。这与多元回归模型中一样， R 2 的增加可以通过增加解 
释变量来实现。相似地，在结构方程模型中，模型拟合的改 
善也可以通过放松参数假定来实现。 

修正指标 

如果模型的拟合不好，可以通过增加额外的参数限定来 
改善模型。一种确定所有增加的参数的方法是基于每个参 
数与拟合函数的差值来确定的 （ S 5 rbom ， 1975)。这种差异 
表明如果增加了这个参数，卡方将会如何减少。这种方法的 
局限在于，如果增加了差异较大的参数,它的估计值可能很接 
近要估计的参数。如果是这种情况，整体模型拟合的改善可 
能很小。约瑞斯科和索尔波姆 ( J 6 reskog & Sorbom , 1981) 
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提供了解决这种问题的模型修正指标。给定参数的修正指 
标等于增加该参数所减少的卡方值。他们建议每次只增加 
一个参数，因为增加一个参数可能减少或消除增加另一个参 
数所带来的拟合度的改善。每步所要增加的参数需要具有 
最大的修正指标.而且应该具有较大的实际意义（更多讨论 
参见《验证性因子分析》)。 
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协方差结构 模型: USREL 导论 


第8节 I 结构系数的解释 

IRMIIIIiaHliaHliaaBMaMRMaMHHHflHNMMKHKaRWHBBMnMflmHI 


B 和 r 中的结构系数可以解释为对内生变量的直接影 
响。 ft 表明，在控制其他变量的情况下，内生变量—个单 
位的变化会导致 7 , ft 个单位的变化。类似地，表明在控制 
其他变量的情况下，外生变量？ j 一个单位的变化会带来内生 
变量卩义,个单位的变化。如果变量已经标准化，解释也应作 
岀相应的调整，在其控制其他变量的情况下，内生变量7 — 
个标准差单位的变化会导致7#,个标准差单位的变化;在控 
制其他变量的情况下，外生变量$ —个标准差单位的变化会 
带来内生变量 7i yi, 个标准差单位的变化。 

这些解释都假定其他变量不变。然而，在实际中，一个 
外生变量的变化可能导致其他变量的变化。如例 3.3 .在其 
他变量不变的情况下， y 13 表明研究对象一个单位的 SES ( f 3 ) 
的变化会直接影响其期望（％ )y 13 个单位的变化。但是在这 
个例子中，模型假定$3的变化会导致&的变化（直接通过 
y 23 )，w 的变化会引起％的变化(直接 通过戽 2 )。因此 ，匕的 
变化既会对％产生直接影响，也会产生间接影响。 因此^ 
对71的总效应必须包含直接效应和间接效应。这个例子表 
明，在解释结构系数时，有必要区分总效应、直接效应和间接 
效应。[ 7] 
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在模型 ry G G 中，内生变量之间的直接效应包 

含在矩阵 B 中，外生变量对内生变量的直接效应包含在 r 矩 
阵中。外生变量对内生变量的总效应可以通过方程 = 一 
— Br \ 获得。如果定义 n =( i — Bn 则〜就 
表示外生变量 e 对内生变量彳的总效应。而内生变量之间的 
总效应要更为复杂一些，因为必须要考虑到内生变量之间的 
相互影响。在例 3.3 中，7^通过 外1 影响7/2，反过来卩 2 又通过 
/? i 2 影响彳2。所以内生变量之间的总效应通过 （1 一 B 广 1 一1 来 
获得，而间接效应可以通过总效应减去直接效应获得。关于 
这方面的详细讨论，可以参考格拉夫和施密特 （Gmff & 
Schmidt , 1982)、福克斯 （ Fox ，1980) ，以及约瑞斯科和索尔 
波姆 (Jdreskog & Sorbom , 1981) 的论述。 

下面我们通过例 3.2 和例 3.3 来展示如何估计和检验结 
构方程中的假设，我们鼓励读者运用附录中的数据来再现这 
些例子。 


例 3.3: 参数估计和假设检验 

例 3.3 中的模型 M a S 传统的联立方程系统，可以通过任 
何一种方法进行估计。表 3.2 列出了二阶段最小二乘法 
(2 SLS ) 、有限信息最大似然法 ( LIML ) 、完全信息三阶段最小 
二乘法 (3 SLS )、 完全信息最大似然估计 （ FIML ) 等方法的参 
数估计。前三种估计方法是通过 SAS 中的 SYSREG 程序来 
实现的， FIML 的估计是通过 LISREL 5.0 实现的。从表 3.2 
中可以看出，不同的估计方法所得出的系数是很相似的。 

接下来我们来介绍表 3.2 中通过 F 1 ML 法所估计的参 
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协方差结构 樓型: USRF 1 导论 


数。研究对象的 SES (匕）对其教育期望（ 71 )的直接效应最 
强。在其他变量不变的情况下， ^一个 标准差的变化会引起 
7 i 0.277 个标准差的变化，而且具有统计显著性。如前文所 
述，直接效应并不能代表& 对 71 的总效应，因为&会影响 
%， 而恥又 会影响因此& 对％ 的总效应要大一些，为 
0.296。&和&对 7 ,也有显著的影响，但是要比&的影响弱 
一点 。&和 ％对％的直接效应都不显著。表 3.3 的底部报告 
了模型从的拟合度，自由度为2,卡方值为 1.88( P = 0.39)。 

虽然模型 M a 的拟合度是比较不错的，但我们也有必要 
考虑该模型的拟合度是否可以改善。例如，在的方程中. 
如果释放参数 y 15 和 y 16 中的任何一个，模型还是可以识别 
的，相似地，也可以释放参数 y 12 和 y 21 中的任何一个。需要 
注意的是，如果对模型的改进并无重要的意义，请不要增加 
这些变量。类似地，如果有理论上的依据，也可以限定某些 
参数为 0( 例如 y M = o )。 

这个模型的对称特性提供了另外一种路径，即每个影响 
研究对象期望的变量都会影响到其朋友的期望，没有理由假 
定影响研究对象期望的变量不会影响其朋友的期望。例如， 
研究对象的期望对其朋友的期望的影响与其朋友期望对他 
的期望的影响应该是相同的。这意味着限定如=成 2 (参见 
Joreskog & Sorbom , 1981)。包含这一系数等同限定的模型 
在表 3.3 的模型 M b 中，增加这一限定后，模型增加了 1个自 
由度，卡方增加了 0.24 o /? 21 =/? 12 的假设可以通过比较模型 
M a * M b & 卡方差异进行 检验; 两个模型之间的卡方差异为 
0.24,自由度差异为1，我们不能拒绝私 1= i 9 12 的假设。因此 
增加如=^12的限定或许是合理的。 



第 3 章结构方程模型 


55 


表 3.2 同辈影响期望结构方程模型不同估计方法的比较 


系数- 

2SLS 

LIML 

3SLS 

FIML 

效应 

t 值 

效应 

t 值 

效应 

t 值 

效应 

z 值 

(3\2 

0.158 

1.39 

0.158 

1.37 

0.157 

1.38 

0.157 

1.38 

y\\ 

0.196 

4.20 

0.196 

4.20 

0.193 

4.22 

0.193 

4.20 

r\2 

0.252 

4.86 

0.252 

4.86 

0.254 

4.97 

0.254 

5.00 


0.277 

5.37 

0.277 

5.36 

0.277 

5.36 

0.277 

5.35 

ru 

0.050 

0.83 

0.050 

0.83 

0.050 

0.83 

0.050 

0.83 

/?21 

0.240 

1.90 

0.240 

1.90 

0.238 

1.89 

0.238 

1.91 

723 

0.047 

0.76 

0.047 

0.75 

0.047 

0.76 

0.047 

0.76 

yiA 

0.254 

5.33 

0.254 

5.33 

0.250 

5.27 

0.250 

5.21 

72S 

0.323 

6.28 

0.323 

6.28 

0.333 

6.53 

0.333 

6.80 

726 

0.201 

4.58 

0.201 

4.58 

0.189 

4.40 

0.189 

4.32 


注:关 于系数的具体意义请参考图 3.1 和文中解释。所有的系数都经 
过标准化 处理。 2 SLS 代表二阶段最小二乘法， LIML 代表有限信息最大似 
然法 ，3 SLS 代表三阶段最小二乘法 . FIML 代表完全信息最大似然法。 
2 SLS 、 LIML 和 3 SLS 是通过 SAS 中的 SYSREC 进行估计的 . F 1 ML 是通过 
LISREL 5.0 估计的。 


表 3.3 同辈对其期望影响的 FIML 估计 


系数 

模型 M , 

模型 Mb 

模型 

效应 

z 值 

效应 

z 值 

效应 

£值 

^\2 

0.157 

1.38 

0.194 8 

2.25 

0.192* 

2.24 

/II 

0.193 

4.20 

0.193 

4.21 

0.192 b 

6.01 

712 

0.254 

5.00 

0.250 

5.00 

0.269 c 

8.13 

/13 

0.277 

5.35 

0.270 

5.42 

0.264 d 

7.40 

7u 

0.050 

0.83 

0.038 

0.69 

0.052* 

1.20 

^2\ 

0.239 

1.91 

0.194 - 

2.25 

0.192 a 

2.24 

/23 

0.047 

0.76 

0.062 - 

1.19 

0.052^ 

1.20 

yi\ 

0.251 

5.21 

0.256 

5.49 

0.264 d 

7.40 

725 

0.332 

6.80 

0.337 

6.93 

0.296 c 

8.13 

726 

0.189 

4.32 

0.188 

4.31 

0.192 b 

6.01 

P 方值 

1.88 

2.12 

4.28 

df 

2 


3 


7 


P 值 

0.39 

0.55 

0.75 


注:罗 马字母表明限定参数相等。关于系数的具体意义请参考图 3.1 和 
文中解释。所有的系数都经过标准化处理。 FIML 是通过 LISREL 5.0 估 
计的。 



56 


协方差结构 模型: LLSREL 导论 


同样，也可以假定研究对象的 SES 对其期望的影响与其 
朋友的 SES 对其自身的期望的影响是等 同的; 研究对象智力 
水平的影响对其期望的影响与其朋友智力水平对其自身的 
期望的影响也应该是等同的;研究对象父母期望对其期望的 
影响与其朋友的父母期望对其自身期望的影响也应该是相 

等的 。 BP 7ll = 726 * 7\2 = 72S * 7l3 = 724 > 7u = /23 o 表 3.3 中 

的模型列出了参数估计的结果。增加了这些限定后，模 
型的拟合度有了显著的改善，即自由度增加了 4个，卡方增 
加了 2.16,考虑到增加这些限定对模型拟合改善的作用，我 
们可以合理地认为，模型是最好的模型。// 

在这个例子中，我们可以看到应用协方差结构估计软件 
的优势在于检验系数等同的限定。下面我们通过例 3.2 来说 
明允许部分误差项相关的限定。 

例 3.2: 参数估计与假设检验 

表 3.4 包含了关于心理失调的三个模型。前文已经证明 
第一个模型 M a (假定 唯 =/? 51 =0) 是可以识别的。本章 
所得到的结果与惠顿 （ Wheaton , 1978： 394 一 395 ) 原文中更 

为复杂的协方差结构模型的结果是相似的。通过 系数如 = 
0.877 和/? 42 二 0.844 可以看出心理失调的高度稳定性 。 7ll = 
0.294 表明研究对象父亲的社会经济地位对其初始的社会经 
济地位(％)有显著的影响，但是其影响在之后消失了。你 3 = 
0.495 也表明了心理失调的高度稳定性。除了研究对象早期 

)对其早期的心理失调症状有显著的影 
响（％)，研究对象的社会经济地位对心理失调的影响均不显 
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著。需要注意的是，这些系数并没有经过标准化处理，不能 
直接比较系数的大小，关于这些结果的更多讨论，参见惠顿 
的研究 ( Wheaton ， 1978)。 

模型 M a 的卡方为47.9,自由度为2,这些指标说明 M a 
的拟合度是比较差的。模型中的每个参数都有理论上的依 
据，因此不能删除。模型 M b 限定早期社会经济地位对其后 

社会经济地位的影响是相同的（私 1 [8] ，增加这一限定 

后，模型增加了 3个自由度，卡方为49.3。因为模 型以 1 是嵌 
套于 M a 的，所以可以通过卡方差异来检验/? 21 =^ 2 的假设， 
结果表明卡方增加了 1.36,自由度增加了 1，因此不能拒绝 
/?21 = /?42 的假设。 

模型 M a 的卡方为47.9,自由度为2,这些指标说明 M a 
的拟合是比较差的。模型中的每个参数都有理论上的依据， 
因此不能删除。模型 M h 限定早期社会经济地位对其后社会 
经济地位的影响是相同的 （fti U ， 增加这一限定后，模型 
增加了 3个自由度，卡方为49.3。因为模型 M h 是嵌套于 
的，所以可以通过卡方差异来检验外, 味 2 的假设，结果表明 
卡方增加了 1.36,自由度增加了 1，因此不能拒绝如=仏 2 的 
假设。 

既然没有看起来合理的其他限定，我们或许应该考虑增 
加参数。模型，修正指标显示，如果模型是可以识别的，释 
放少 24 可以显著改善拟合度。增加参数少 24 的结果在表 3.4 中 
的模型叭中。释放了 ％,后，模型依旧是可识别的，并且拟 
合度有了显著改善。少 24 是预测早期的社会经济地位（％)和 
预测晚期的社会经济地位（％)的方程中误差项的协方差。这 
种误差项的相关也被称做序列相关 （Kessler 8^ Greenbrerg , 
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1981)，它表明模型中可能没有包含同时影响两个因变量的 
其他自变量。例如，布劳 （Blau) 和邓肯的研究表明教育获得 
会影响到社会经济地位，而在预测&和 73 的方程中都没有 
包含教育获得变量，这会导致^和 G 的相关。模型 H. 的 
卡方值为1.58,自由度为2,这表明模型 H, 的拟合度是非常 
不错的。// 


表 3.4 心理失调结构方程模型的 FIML 估计 



模型 

模型 M h 

模型吣 

乐奴 

效应 

z 值 

效应 

z 值 

效应 

z 值 

7 \\ 

0.294 

6.52 

0.294 

6.52 

0.294 

6.52 

卢21 

0.877 

44.10 

0 . 861 3 

59.11 

0 . 896 * 

69.45 

721 

0.026 

1.15 

0.031 

1.38 

0.021 

0.92 

/?31 

- 0.010 

3.80 

- 0.010 

3.90 

- 0.010 

3.72 

ysi 

— 0.002 

0.76 

- 0.002 

0.73 

_ 0.002 

0.78 

/?42 

0.844 

40.01 

0 . 861 - 

59.11 

0 . 896 " 

64.95 

夕43 

— 0.057 

0.18 

_ 0.035 

0.11 

0.018 

0.06 

741 

0.060 

2.50 

0.055 

2.33 

0.046 

1.93 

/?52 

— 0.002 

0.73 

- 0.001 

0.58 

- 0.002 

0.88 

卢53 

0.495 

14.95 

0.495 

14.96 

0.495 

14.92 

• 

ysi 

- 0.003 

1.34 

— 0.003 

1.38 

一 0.003 

1.31 

023 

1.69 

2.68 

1.69 

2.68 

1.72 

2.72 

045 

2.37 

4.28 

2.37 

4.28 

2.35 

4.41 

024 

— 

—— 

—— 

— 

一 35.31 

6.1 

卡方值 

47.95 

49.31 


1.58 


df 


2 

3 


2 


P 值 

0.000 

0.000 

0.454 


注:罗 马字母表明限定系数相等。对于这些系数的具体意义请参看图 
3.1 或者文中关于这些系数的解释。所有的系数都没有标准化，这些估计采 
用了 LISREL 5.0 的完全信息的最大似然方法。 
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第9节 | 小结 


在协方差结构模型的框架下，本章所介绍的结构方程模 
型是很有意义的。第一，协方差结构估计软件对于估计那些 
非标准化的不含测量模型的结构方程模型是很有帮 助的; 第 
二，应用协方差结构模型时，必须在潜变量之间构建结构方 
程模型，这与本章所讨论的结构方程模型是一样的，不同的 
是在协方差结构模型中，结构方程模型是在潜变量之间建立 
的，而本章是在观测变量之间构建的。我们将在下一章中讨 
论如何将第2章中的测量模型和本章中的结构方程模型整 
合在一起。 




mM 
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因子模型在没有考虑潜变量之间关系的情况下从观测 
变量中估计潜变量，可是潜变量之间的关系通常是研究者的 
兴趣 所在; 而结构方程模型虽然关注变量间的结构关系，却 
必须假定所有的变量都不含测量误差，而因子模型技术表明 
这种假设往往是不现实的。整合了因子分析和结构方程分 
析方法的协方差结构模型在一个模型中同时从观测变量中 
估计潜变量，以及在潜变量中估计结构关系。因此，它既克 
服了两者的弱点，又吸收了两者的优点。既然我们已经介绍 
了测量模型和结构方程模型，再介绍协方差结构模型就相对 
简单了。在本章中，读者需要了解那些连接这两个部分的假 
定以及如何将我们前面关于模型的识别、参数估计及假设检 
验应用到协方差结构模型中。 


第 1 节 I 数学模型 


第4瘫协方差结构楱型 


协方差模型的结构部分与第3章中讨论的结构方程模 
型的数学表达是一 致的： 

7|= B?y + r 4 + C [4.1] 

t ] 为 (rX 1) 维的内生潜变量的向量集 j 为 （sX 1) 维的外生 
潜变量向量集，《为 （rX 1) 维的方程误差的向量集。 B 为 
( rXr ) 维的内生变量之间关系的系数矩阵 ， r 为 （rX s ) 维的 
外生变量与内生变量之间关系的系数矩阵。如果将6定义 
为 （1 — B ) ，则方程 4.1 可以写做 Bri = r 4 +《。 

第3章中关于结构方程模型的假定在这里依然适用。 
第一，假定所有变量的均值为0,即 E ( ry ) = E (《） = 0, 
E (^)= 0；第二•没有冗余的方程，即 ( I - B )- 1 = fi 1 是存 
在的； 第三，方程的误差项与内生变量不相关，即 E (^) = 
0或= 0。表 4.1 列岀了与协方差结构模型的相关 
假定。 

关于变量间协方差的定义在这里也适用。外生变量的 
协方差矩阵定义为. •4> = E ( K )。 方程的误差项矩阵定义为 
平为对称矩阵，但不一定是对角线矩阵。由方程 
4.1 及关于模型的相关假定可以得出内生变量的协方差矩阵 
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CC ) V ( ry ) = E ( 7 , 7 / ^) =B 1 ( + (推导过程可以参见 

第3章)。 


表 4.1 协方差结构模型 


矩阵 

维度 

均值 


协方差 

维度 

注 释 

V 

(rX 1) 

0 

COV ( i /) = ECt / tj ) ( rXr ) 

内生潜变量 


(sX 1) 

0 

中= 

E (^) 

(s X s ) 

外生潜变量 

( 

(rX 1) 

0 

平= 


(rX r ) 

方程误差 

B 

( rXr ) 

— 


—— 

—— 

内生变量间的直接效应 

•• 

B 

(rX r ) 

— 


— 

-— 

(卜 B ) 

r 

f r V c 、 





外生变量对内生变量 

1 

v r /\ S / 





的直接效应 

X 

(qX 1) 

0 

Zr , 

= E ( xx ’） 

(q X q ) 

外生观测变量 

A . 

(q X s ) 

— 


— 

— 

• T 在^上的因子负载 

S 

(qX 1) 

0 

0^ = 

= E ( 8 S ，) 

(q X q ) 

X 的独特因子 

V 

( rXl ) 

0 

CQV ( iy )= E ( i 7 rj ， ) ( rXr ) 

内生共同子 

y 

(pX 1) 

0 

hyy z 

= E ( yy ') 

(p X p ) 

内生观测变量 

Av 

(pX r ) 

— 


— 

_ 

: y 在上的因子负载 

€ 

(pX 1) 

0 

0,= 

E (€€) 

(pX p ) 

y 的独特因子 


注 :结构方程: 
因子 方程： 

协方差方程： 


t^Btj 十 K+f 
X = iVrg +6 

y = A v rj-fc 



"Avfi-u ror' + 平 ) 6' - 1 a: v +0 

AvB- j rOA ； ( 1 


_ A , Or B ^Ay 

\u OAj - ! - _ 


[4.1] 

[4.2] 

[4.3] 


[4.4] 


假设： 

(1) 所有变量的均值为 0: E(rj) = E(f) =0； E(§) =0, E(x) = E(8) 
= 0 ； E(y) = E(f) = 0 ： E(-q) — 0 ； 

(2) 共同因子与独特因子不 相关： E (鉍 '）= 0 或 E(64 f ) = 0 ; E(tk' 〉= 
0 或 E(oj ’） = 0; 

(3) 不同方程中的独特因子不 相关： E(d') = 0或 E (€ S ) = 0； 

(4) 方程中的内生变撤和误差项不 相关： E(G') = 0或 E(g') = 0; 

(5) 没有冗余的方程： (1 —B 广 1 存在。 
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与第3章中的结构方程模型不同，协方差结构模型中的 
T ? 和4都是不可观测的潜变量。 ry 和$通过一对因子模型与 
观测变量 x 和 y 连接起来。 


A . j + 6 

[4.2] 

A v tj- he 

[4.3] 


x 是 (qX 1) 维的外生观测变量集， y 为 （pX 1) 维的内生观测 
变量集， Ai 是 ( qXs ) 维的观测变量 j : 在潜变量$上的因子 
负载矩阵， A v S ( pXr ) 维的观测变量 j ； 在潜变量 fy 上的因子 
负载矩阵和 f 分别为 (qX 1) 维和 （pX 1) 维的独特因子向 
量集。在每个因子模型中，独特因子之间是可以相关的，也 
就是说 COV (5) = E (55') =0 5 和 COV (€) = E(€€)=e 是对 
阵矩阵，但不一定为对角线矩阵。其实方程 4.2 和方程 4.3 
与第2章中的方程 2.1 和方程 2.2 是一样的。 

第2章中关于测量模型的假定依然适用于协方差结构 
分析中的测量部分。第一，同一方程中的共同因子与独特因 
子不相关，即 E ( gS ，) = E ( tj € )=0； 第二，不同方程中的共同 
因子与独特因子不相关，即 E (^ ， )= E ( tjS ， )=0； 第三，误差 
项占，（，？之间不相关，即 E ( Se f )= E ( S ^)= E ( e ； ， )=0 o 
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第 2 节 I 协方差结构 

mi 麵 _|— wii 麵 _1 


既然假定所有变量的均值为0,观测变量的协方差矩阵 
就可以定 义为： 




■ fmm mb f — 

y_ 

■ X 」 L X _ 


yy yx 


是 （p + qX 1) 维的向量•将方程 4.2 和方程 4.3 代入， 


得到: 


I： = E 


"(A v ty 

_( A ^ 


+ f)(A v ?y+f)’ I (A v ?y 4~£)(八』《 +6) ’ 
十 5)(A v i|+£)’ (A_r 《+5)( A.:4+5): 


通过乘法得到: 


" A . vT 7? y ， Av +^ / 
+ \ yTje r + ft]'\ f y 

A ., gTy ' A' y + 5 〆 
_+ Aj §£’ + St] f \ f y 


A. v T|4A ， r+€8 / 

+ A v tj 6 +4’ A : 

+ A^5 f + 5$X. 


通过数学期望的变换及协方差为0的假定，可得如下 
方程： 

fA v B (ror + ^)B ，] A； + 0. A V B 1 roA.l 


A,Or'6’ _1 A’ v 


A r OA x +0, 


[4.4] 


第 4 素协方差结构模型 


方程 4.4 通过模型参数来表达观测变量的方差和协方差。假 
定模型是可识别的，估计的过程就是根据方程4. 4 找到八个 
参数矩阵的值来产生一个与样本协方差矩阵 S 最为接近的 
协方差矩阵2：。 
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第3节 I 协方差结构模型的特例 




将验证性因子分析模型和结构方程模型作为协方差结 
构模型的特殊情况是非常有意义的，因为这样不仅可以更好 
地理解协方差结构模型，还可以理解如何运用协方差结构模 
型的软件来估计验证性因子分析模型和结构方程模型。 


因子分析模型 

通过限制 B = o ， r = o ， 平 = o ， a v = o , 0 = 0 ,协方差结 
构模型的结构部分简化为《 =0(( 的方差也为 0) ，^变量的 
测量部分简化为这样 I 变量的测量部分与就与方 
程 2.1 等同： x =\^+ S 0 控制这些限定条件后，方程 4.4 就 
简化为 = 这与方程 2.3 中的左上部分是相 

同的。 

相似地，通过限制 b=(k6 = i )， r = o , a v = 0 , 

0 =0,以及把 I 变量的测量部分删除，结构方程模型就变 
成了 r / =《， n 的协方差为平。此处的平就相当于 T 变量因 
子模型中的 中， 控制这些限定条件后，方程 4.4 就简化为2：= 
A v WA' v +0 。 
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结构方程模型 


通过限定潜变量与观测变量相同，就可以将协方差结构 
模型中的测量部分删除，从而得到结构模型。限定 A . =1， 
a =0, A f =1，0 =0之后，方程 4.1 和方程 4.2 就简化为\= 
^ y = fi ， 即潜变量与观测变量是相同的，方程 3.1 和方程 4.1 
都保持不变。限定& = 1，03 =0, =1，0, =0之后，就可 

以通过方程 4.4 得到方程3.5。因此，结构方程模型就是协方 
差结构模型的特例，在施加参数限定后，就可以运用协方差 
结构的软件来估计结构方程模型。 
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第4节 I 实例 


下面我们通过两个同时包含了测量部分和结构部分的 
例子来说明模型的设定、识别、估计及检验^ 

例 4. 1:含有测量误差的追踪数据模型 

本例将《验证性因子分析》中的例1整合到了第3章例 
3.2 的追踪数据模型中，尽管有些简化，但基本上和惠顿 
( Wheaton , 1978) 的原模型是一致的。如图4.1，该模型包含 
一个外生变量和五个内生变量。父亲的社会经济地位（&)、 
研究对象三个时期的社会经济地位（％，％和&)，以及研究 
对象两个时点上的心理失调症状（％和 rj ,) 0 结构模型与例 
3.2 的最后一个模型 一样： 


V 


"0 

0 

0 

0 

0~ 




>]|" 





夕21 

0 

0 
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首先假定方程的误差项只在同一时点上是相关的，即在时点 
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2上预测研究对象社会经济地位和心理失调状况的两个方程 
的误差项是相关的，在时点3上，预测各个变量的方程的误 
差项也是相关的。因此少矩 阵为： 



时点 2 〜 6 °' 7 时点 3 


注: FSES 为父亲的社会经济 地位; SES 1 为研究对象时点1时的社会经 
济 地位; SES 2 为研究对象时点2时的社会经济 地位; SES 3 为研究对象时点 
3时的社会经济 地位; PSY 2 为时点2时心理失调症状的数目 t PHY 3 为时点 
3时心理失调症状的数目； PHY 2 为时点2时精神生理失调症状的数 B ; 
PHY 3 为时点3时精神生理失调症状的数目。 

资料 来源： Wheaton , 1978。 

图 4.1 心理失调症状的协方差结构模型 


与例 3.2 不同的是，该模型包含了测量模型。外生变量 
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父亲的社会经济地位 (&) 是没有测量误差的，即 A =匕。 
A ., = I ，队= 0。内生变量的测量误差比较复杂，研究对象 
三个时期的社会经济地位（^， 72 和％)是没有测量误差 
的 ， BP % ， V 2 = 74 = y -^ 也就是说 Afi = = 

= 1，〜4 = I = 0。内生潜变量心理失调是通过 
心理失调症状的数目和精神生理失调症状的数目两个观测 
变量来测量的，即: y 3 和: yi 在潜变量％上的因子负载为 
和，>和 h 在潜变量彳 5 上的因子负载系数为 A 〖 5 ，和 
A &。 分别设定因子负载 A 3 V 3 和 A 〖 5 的值为 1( 参考验证性因子 
分析第2章）。由于 COVCi ^ S 模型中最重要的参数——因 
为 COV ( 17 ) 要通过其他参数矩阵来表达，因此矩阵内生变量 
7只能通过限定因子负载矩阵 A v 来设定，所以^的测量模 
型可以记做： 



与《验证性因子分析》中的例1相同，模型假定同一变量的测 
量误差在不同时点上是相关的 ，即& 关 0 , 參 0 。 由于 

: Vl ， h 和: y5 不含测量误差，因此 =022 =^55 =0, 0 , 的协 
方差矩阵可 写为： 
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因此图 4.1 中测量模型有两种表达方式 :第一 ，当没有测量误 
差时，观测变量和潜变量是相同的，都包含在椭圆 形中； 第 
二，当包含测量误差时，观测变量在矩形中，而潜变量则在橢 
圆形中，通过单向箭头将观测变量和潜变量连接起来。// 

例 4.2: 包含测量误差的非递归模型 

本例（图 4.2) 在例 3.3 中结构模型的基础上增加了内生 
变量的测量模型。 

在这个例子中，外生变量不含测量误差，即 x = ^ A , = 
1 ，蛛=0,而每个内生变量都有两个测量指标。内生变量 71 
(研究对象的期望)是通过力（研究对象的教育期望)和研 
究对象的职业期望)来测量的。同样，内生变量％ ( 研究对象 
朋友的期望)是通过％(研究对象朋友的教育期望) 和％ (研究 
对象朋友的职业期望)来测量的。因此内生变量的测量模 型为： 
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协方差结构 模型: USREL 导论 



注: R - PASP 为研究对象父母的 期望; R - IQ 为研究对象的 IQ ; R - SES 为 
研究对象的社会经济 地位; F - SES 为研究对象朋友的社会经济 地位; F - IQ 为 
研究对象朋友的 IQ ; F - PASP 为研究对象朋友父母的 期望; R - EASP 为研究 
对象的 期望; F - ASP 为研究对象朋友的 期望； R - OASP 为研究对象的职业期 
望; F - OASP 为研究对象朋友的职业期望。 

资料来源 ： Duncan et al . ，1971。 

图 4.2 同辈对期望影响的协方差结构模型 

模型设定因子负载 A 和 A 〖 2 的值为1。邓肯等人 （ Dimcan ， 
1971) 假定测量误差之间存在如下 关联： 



‘和0 24 (对应地〜和化 2 )表明研究对象及其朋友的职业期 
望的测量误差是相 关的; 研究对象及其朋友的教育期望也是 
相关的。由于类似的测量误差可能同时存在于对研究对象 
和其朋友职业期望的测量中 ，上 述误差相关的设定是合理 
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的。〜（等于％ )和心（等 于心） 表明研究对象的教育期望 
和其朋友的职业期望是相关的（反之亦然），这看似是没必要 
的；如果教育期望和职业期望的测量误差是相关的，假定其 
在同一人身上是相关的则是很合理的（比如〜和•因此 
在我们的分析中， ^和礼 被限定为0。// 
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协方差结构 横型: LISREL 导论 


第5节 I 模型的识别 


与因子分析和结构方程模型一样，在估计协方差结构模 
型之前，必须要确定模型是否可以识别。协方差结构模型的 
识别与我们之前的讨论并没有什么不同，因为其增加的只是 
参数的个数。 

模型可识别的一个必要条件就是所有估计的自由参数 
的个数必须小于或等于样本观测变量协方差矩阵 S 中非冗 
余元素的个数。当计算自由参数的个数时，必须考虑到参数 
等同限定及对称矩阵等因素。例如，限定 yu = y 12 = y 13 =yM 
后，自由参数的个数为1，这是因为知道了其中的一个参数 
后，就可以计算出其他参数。再比如，少 12 和 少 21 只能算做一 
个自由参数，因为在对称矩阵中，平 12 =平 21 。如果用 t 来表 
示协方差结构模型中自由参数的数量，则模型可识别的必要 
条件为： t < 34(p + q )(p + q 十 1) 。 

尽管前文讨论了判断特定模型是否可识别的条件，然而 
这些条件并不能简单地用来判断饱和协方差结构模型是否 
可以识别。通常，模型是否可以识别必须通过运用观测变量 
的方差和协方差来证明模型中的参数是可以求解的。对协 
方差结构模型而言，首先要证明测量模型中的参数是可以识 
别的，然后运用因子的方差和协方差来证明结构模型是可以 
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识别的，需要注意的是，协方差分析中的结构方程模型的识 
别是不能直接运用观测变量的方差和协方差来证明的。因 
子之间的协方差可以通过观测变量的方差和协方差来解决, 
这是证明测量模型可以识别的充分条件。需要注意的是，除 
非首先证明测量模型是可以识别的，否则运用因子间的方差 
和协方差来解决结构参数并不是证明结构模型可以识别的 
充分条件。 

验证一个模型是否可以识别是一件非常困难的事情。 
约瑞斯科和索尔波姆认为可以通过计算机估计软件来证明 
一个模型是否可以识别，因为在对参数进行最大似然估计 
时，估计软件会计算信息矩阵 （Jdreskog Sorbom , 1978)。 
简单地说，信息矩阵就等于参数的方差和协方差矩阵，约瑞 
斯科和索尔波姆提到，如果信息矩阵是正定的，那么模型几 
乎都是可以识别的，如果信息矩阵为非正定的，则模型是不 
可识别的。 

读者需要注意的是，这并不是模型可以识别的充分条 
件。如果信息矩阵是正定的，模型也有可能无法识别或者估 
计的参数是没有意义的。另一方面，如果通过观测变量的方 
差和协方差证明模型是可以识别的，信息矩阵也有可能是非 
正定的，或者在运行程序时会给出模型不可识别的错误指 
示。如果要用信息矩阵来判断一个模型是否可以识别（我们 
并不推荐这种方法），研究者应该采用约瑞斯科 （ Mreskog ， 
1978) 的程序。假如根据信息矩阵显示模型是可以识别的， 
读者应该尝试设定不同的初始值，如果得到的估计结果相 
同，研究者会有较大的信心(仍然不确定)认为模型是可识别 
的。可是当条件几乎满足时，如果模型仍然是不可识别，那 
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协方差结构 模型: LLSREL 导论 


么其结果也是没有意义的。本书推荐研究者通过用观测变 
量的方差和协方差矩阵来解决模型参数的方式来证明一个 
模型是否可以识别。 

例 4.1 的识别 

这个例子中的结构模型与例 3.2 是相同的（比较图 3.1 
和图4.1)。不同的是例 3.2 中的结构方程的变量是可观测变 
量，而本例中有些变量是潜变量。在例 3.2 中，我们通过每个 
参数(/3, y ， 氺)都可以运用观测变量的方差和协方差（卩，幻 
来解决的方式证明模型是可以识别的。既然例 3.2 与本例中 
的结构方程模型是相同的，其中关于每个参数的代数表达就 
可以通过潜变量的方差和协方差来解决，如果潜变量的方差 
和协方差是可以识别的，那么本例中的结构模型也是可以识 
别的。 

要识别潜变量的方差和协方差，就必须证明其中的每个 
参数都可以用观测变量的方差和协方差来解决。由于观测 
变量和潜变量是通过测量模型联系在一起的，我们就必须从 
测量模型开始。 

^ 1 — y\ = y]\ yi — yji 

yz = T]3 ^3 3^4 = A 43 7^3 + (4 3^5 = T]\ 

+ f 6 yi = ^75 TJr } + c 7 

我们可以通过乘以每对测量方程，取期望，然后运用某些变 
量不相关的假定（方程中用线划去的部分）得到观测变量 X 
和3^的协方差。观测变量的方差为： 
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~ ^11 

Cfylyl = VARC^l ) 

a y2 y2 = VAR (7；2) 

Oyiyi = VAR( 7^3 ) + ^33 + 2COV (”” (T^ 

Gyiyi =A4 3 " VAR(^ 3 ) + d\ A + 2A;J 3 Covers, d ) 

(3y5y5 =VAR(^ 4 ) 

0# = VAR( 75 ) + 汐 ‘66 + 2COV(t75 i~r^ 
o y 7y7 =A 7 v 5 2 VAR(y) +^77 + 2A 7 S Covers v 
观测变量的协方差为： 

CTjlyl = COV(f 1 9 7j]) 

(jjj y 2 = COV(f 1 » 7 ^ 2 ) 

(7x1^3 = COVC^! ， 7/3) + COV($l ， r 3 -) 

叫 4 =Aj 3 COV(e l , ?3 )+ COV( 。， rr^ 

Os\yS = COV(fi ^ Tj.\ ) 

(j rlv6 =c()v(f!, 75 ) +COV(Ci 

<1x1^7 =A 7 5 COV( 6 !, 75 ) + COV(f), rr^ 

^2 = COV( 771 , 72 ) 

(Jvi.v3 = COV (rj' ， 7 / 3 ) + COV( 刃， h) 
aw =A4 3 COV(7 ；i , 7 / 3 ) + C()V(, r 4 > 

0y\y5 =COV(7Ji ， 7/4) 

o y \ye =COV(^i, 7 / 5 ) + COV(^i ^ nr) 

0y\ y7 =A 7 5 COV(^i , 7j- 0 ) + COV( J 71 , n ) 



协方差结构 横型: LISREL 导论 


Oy2y3 ~ C()V( Tjz » ”3 ) + COV( 7]2 i T3 ) 

~ A J 3 CC)V( 7^2 ♦ 7^3 ) + COV< 72 、 「i ) 

Oy2y5 = COV( )^ 2 , Yj\) 

Oy2y6 =COV(^2 » ^5 ) + COV( 7 ^ Tg^ 

av2>7=A? 5 COV( 72 , ?5 >+COV< 72 v h) 

a y zy\ = A43 VAR(w) + @^ + COV( % ， 14) 4 * A43COV( » ( 3) 

O y 3 y 5 = COV( 7 3, w)+ COV < T4 — tt > 

〜 6 =c()v( 7 ” 7 5 ) + g 36 +COV( 7 ^ Q> + COV( 7 ” 

Oy3y7 ~ A^CC^V C TJ^ » 7^5 ) ~h @57 H - COV( 如， < 7 ) 

+ A《 5 GO¥Vt^4 
= A ; 3 COV (7^ , 7^ 4 ) + COV (^ , r .i ) 

^ v 43-6 ~ A^COVC 7 ^ 3 , ^5 ) "I - A 43 COV( Tjz ^ «6 ) 

+ GOVirf ^ rr ^ 

(jy^y? = A "43 A 75 COVC 7^3 » r/s) + ffp + A 43 COV( 73 ， i 7 ) 

+ A? 5 COV( 75 rfr> 

^y5y6 =cov( 7 ” 73 )+ cov( ? 4 ，“) 

a y5 y7 =A《 5 COV(々4 ， r /^) + COV< tj- a —rr) 

W =A^ 5 VAR( rj 5 )+^+ COV( 7 ^ o) + A? 5 COV( 7 s ，n> 

这些方程看起来很复杂，我们可以通过误差项与因子不 
相关及某些误差项不相关的假定来简化这些方程。运用这些 
假定后，我们很容易发现许多参数是可以识别的，因为它们 
都可以用观测变量的协方差来解决。例如 var ( 7i )， 
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var ( 72 )， var ( 74 )， 以及除 cov ( 73 , ％)之外的潜变量的 
协方差都是可以识别的。 

因子负载系数可以通过下列方程证明是可以识 别的： 

A y • 

43 = ^x\yi ~ < 7 xl >3 

A 75 = CT x l y 7 ~ (Jily6 

由于因子负载系数是可以识别的，其余的潜变量之间的协方 
差也是可以识别的： 

COV ( 7^3 » r / s ) — o y 3 y 7 +入 ;5 

剩余的参数可以通过下列方程证明其是可以识别的： 

VAR(? ； 3 ) +又匕 

036 = ^y3y€ — COV( "3 ，％ ) 

^47 = ^ y 4 y 7 — 又又 ? 5 COV ( %， ) 

V AR ( 7^5 ) = CTy 6 y 7 A 75 

^33 ~ (Jy 2 y 3 — VAR( 7/3 ) 

0 \i = Oy\y\ — A : 3 2 VAR( 恥 ） 

^66 = a >6^6 _ VAR( 7^5 ) 

Q'rj — O ylyl — A 75 2 VAR( 7 /s ) 

测量模型中的所有参数以及潜变量的方差和协方差都已被 
证明是可以识别的。因此可以通过例 3 . 2 中的方式来证明结 
构方程模型也是可以识别的。// 

例 4.2 的识别 


如上文所述，如果潜变量的方差和协方差是可以识别 
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协方差结构 横型: USREL 导论 


的，那么该模型的结构部分也是可以识别的。本例的证明过 
程与例 3.3 的证明过程是一样的。因为 x = 矩阵就等于 
•T 的方差和协方差， 因此中 是可以识别的。^的测量方程 
如下： 


3^1 = r j \ - I - ^1 = ^' l \ 71 

yz — t ]2 (z yi 7 2 

因为 f 与< 不相关，因此很容易证明 7 和 f 的协方差是可 
以识别的。例如，将 y , 的方程乘以匕.然后取期望得到 
E ( y ] X ]) =E(^iXi) + E((] j：]) ,或者 ( j y \ T \ =COV( t 7 i ，： r!) = 
COV( fl ，&)。 ： y 的协方差方 程为： 


ay ^ y \ = VAR ( j ^ ) +6 ' n 

^ y 3>3 ~ VAR ( rj 2 ) ^33 

vl . v 2 = A'21 VAR ( 7 ]\ ) 

[ 2 C ()V(〜 ， * ) 


a y 2y2 =A2rVAR(77!) +d\ z 
<T.v4v 4 =A f VAR(7； 2 ) +〜 

<7.vb-3 = COV ( 7Ji t 7 / 2 ) d\ 3 

fy y 2\3 = A ^ C ( ) V ( 7j” r/2) 


o y 2yA = A 21 A 7^1 J 72 )+ 0*24 = 义匕 VAR( Tjz ) 


方程中已经删除了假定为 0 的协方差。由于10个方程中含 
有11个自由参数，因此参数是不可识别的。如果将^的测 
量模型与其他模型单独分析•模型也是不可识别的。但是在 
协方差结构模型中，每个参数的识别都会用到整个模型的信 
息。在这个例子中，:^，3. 4 , x, 的方程可以证明7的测 
量模型是可以识别的。其协方差方 程为： 


<7 必 1 = 又 L COV( 】，$〗 ） 

^>4x1 = A:2COV( "2 ， ) 

既然 cov( 7l ， 心）和 COV( 72 , 心） 是可识别的，那么两个因 
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子负载系数也是可以识别的。由于因子负载系数是可以识 
别的，内生变量的方差也是可以识 别的： 

VAR( 7^1 ) = (Jy\y2 A 21 
WAR(rj 2 ) =(T y 3 yA 

这些表明 A , 、 、仏 3 、〜都是可以识别的，下列方程表明其 
余的参数也是可以识 别的： 

COV( 7^1 ♦ 7 2 ^ ~^.vl>4 A 42 

0; 3 = (j v i v 3 — COV (r/\ ^ rjz) 

024 ~ a y 2 yi —^21 A^COVC 771 » TJz) 

每个参数都可以通过观测变量的方差和协方差来求解，因此 
模型是可以识别的。其中有些参数过度识别了，读者应该尝 
试用不同的方法来解决这些参数。// 

证明参数是否可以识别包含了许多步骤、许多参数，许 
多地方都可能会出错。为了避免这些错误，本书建议首先将 
测量模型和结构模型的方程用小写的罗马字母代替大写的 
希腊参数，这样有助于减少出错的几率;第二步•用罗马字母 
列出所有的协方差 方程; 第三步，列出所有需要识别的参数; 
第四步，将这些列表复制 几份; 第五步，当证明了某个参数可 
以识别后，就在列表上划去该参数以及包含该参数的协方差 
方程，直到证明所有的参数都是可以识别的；第六步，稍等一 
天之后，重复上面的五个步骤。尽管这看起来需要大量的工 
夫，但是请记住，如果模型是不可识别的，所有的分析就都是 
没有意义的。 
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第6节 I 模型的估计 

兔8^泌匕碑381浓離 ^ ss's n 版猶 


证明模型可识别之后，就可以估计模型了。协方差结构 
模型可以通过我们在第3章中讨论的完全信息技术来 估计: 
未加权的最小二乘法 ( ULS )、 一 般最小二乘法 ( GLS ) 以及最 
大似然法 ( ML )。 每个参数的估计值都是使得预测的协方差 
矩阵£与样本中观测的协方差矩阵 S 之间的差异最小，这两 
个矩阵之间的差异的定义是由估计方法决定的。与验证性 


因子分析相比，协方差模型的估计并没有复杂多少，不同的 
是估计的协方差矩阵是由下列方程决 定的： 


/\ 



A y B Uror +¥) B’- A J 

八 

-0 

A,B 】roAl 

八 八八 , C , 八！ 

A,Or' B — 

八 /\ 八 1 八 

A., OA x + 0,? 


ULS 估计法可以在不假定观测变量的分布时，得岀一致 
性的估计，然而其缺点是估计值依赖变量的测量尺度而且没 
有相关的统计检验技术。如果观测变量是多元正态分布的， 
那么 GLS 和 ML 的估计都具有渐进性的特征。 ML 最接近 
无偏估计，因为其抽样分布的方差更小，而且接近正态分布。 
GLS 估计与 ML 估计基本相同。关于不同估计方法的特点， 
可以参考《验证性因子分析》中的讨论。 

与因子模型或者结构方程模型相比，协方差结构模型通 
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常有更多的参数。因此估计模型时，其耗费的时间也更长。 
每种估计方法都通过迭代程序寻找更好的估计参数。迭代 
的第一步是要提供初始值，提供了较精确的初始值，会大大 
加快协方差模型估计的收敛 。 LISREL 5.0 的提供了模型初 
始值的选择程序 (Jdreskog & Sorbom , 1981:1.32)。 
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协方差结构 楔型: LISREL 导论 


第7节 I 模型拟合度 


与第3章中的结构方程模型及《验证性因子分析》相比， 
评估协方差结构模型拟合度的技术并没有什么不同。这里 
只是回顾一下前文的介绍。 

检验参数估计值 

协方差结构模型中测量部分的参数检验和解释与因子 
分析模型是一样的。同样，其结构部分的参数检验与解释也 
可以参照结构方程模型。负的方差、大于〗的相关系数及不 
合理的参数估计值都表明模型存在问题。可能的原因有错 
误的模型设定、错误的协方差矩阵、不可识别的参数或者错 
误的模型命令等。 

需要特别注意参数的符号。有些软件将 B 记作 “ BETA ”， 
有些则将 fi 记作 “ BETA ”， 这种情况下，需要注意软件釆用的 
是哪种定义。如果对两个矩阵中的参数施加等同限定，那么 
这两个被限制为“等同”的参数的符号就可能是相反的。例 
如，通过 LISREL 5.0 限定 乂等于 ，从而导致凡 =— y _。 
如果不留意这些区别，就很容易岀错。 




第 4 章协方差结构模型 


87 


参数的方差与协方差 

参数的方差可以通过 GLS 和 ML 来估计。在观测变量 
服从正态分布的假定下，所需估计的参数就接近正态分布， 
从而可以通过 z 检验对单个参数进行检验。通过估计值的 
协方差就可以计算其相关系数，如果两个参数的相关系数很 
大，就表明很难区分这两个参数，这些模型识别不足的讨论 
可以参考《验证性因子分析》。 

卡方拟合度检验 

在观测变量服从正态分布的假定下， ML 和 GLS 估计可 
以提供关于设定模型与零模型之间的卡方检验。相对于自 
由度而言，较大的卡方值表明模型并没有很好地拟合数据。 
协方差结构模型中，卡方检验的自由度为 df = l /2 (p + q ) • 
(p + q +1) — t ， t 为模型中自由参数的数量。 

卡方值的差异检验可以用来比较嵌套模型。如果模型 
M , 可以通过限定一个或以上的模型 M 2 中的参数得到，那 
么模型 Mi 就嵌套于模型 M 2 中。如果模型嵌套于模型 
M 2 中，那么模型的卡方值为 X 〗，自由度为€^ 1; 模型吣 
的卡方值为 X ;，自由度为 df 2; 对于大样本而言， X 2 = XI - 
，服从自由度为 df = dL — df 2 的卡方分布。特定自由度 
下较大的卡方值表明，对于模型施加额外的限定应该被 
拒绝。 
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协方差结构 横型: LISREL 导论 


样本规模对假设检验的影响 

Z 检验和卡方检验都会受到样本规模的影响。 [9] 表面看 
来，样本规模对关于模型拟合度卡方检验的影响与其对单个 
参数 Z 检验的影响是相反的。例如，样本规模越大，模型拟合 
度就有可能越差，而单个参数 Z 检验的统计显著性则会提 
高。然而，这些看似相反的影响，其实都是合理的。 

当样本规模增大时，抽样分布的方差会变小，即由于样 
本规模较大，随机抽取的两个样本之间的差异可能会比较 
小。既然抽样分布的方差较小，那么较差的模型拟合度很有 
可能是因为模型假设是错误的（比如，参数并不等于假设的 
值)而非因为抽样误差造成的（即参数不等于假设的值是因 
为抽样有偏）。因此，当样本规模增大时，参数估计值之间较 
小的差异都会变得显著。 

卡方检验也是类似的情形，可以把卡方检验看做检验模 
型预测的协方差与观测变量的协方差之间差异的一种方法。 
当样本规模增大时，预测协方差与观测协方差之间较小的差 
异都会变得显著。同样，估计的参数值与假设的参数值之间 
的细小差异也会变得显著。 

因此，当样本规模很大时，因为模型的拟合度在统计上是 
不可接受的，几乎任何自由度为正的模型都有可能被拒绝。 
这种情况被本特勒和博内特 (Bender 8^ Bonett , 1980) 称做“最 

小错误”模型。最小错误模型是基于卡方检验而拒绝预测协 
方差与观测协方差之间差异很小的模型。他们认为即便可能 
会有更好的模型来解释这些差异，但最初的模型却提供了最 
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为重要的解释。约瑞斯科和索尔波姆 (Jdreskog & Sorbom , 
1978) 也表达过类似的观点，统计问题并不是假设检验的唯 
一标准，只是提供判断模型是否很好地拟合了数据的参考。 

在协方差结构模型中，样本规模与统计检验的关系其实 
很简单。对 ML 和 GLS 估计来说，卡方统计量为 （N — 1)^， 
N 是样本规模，△是依赖于观测变量协方差矩阵而与样本规 
模无关的量。例如，用同一个模型来估计两个不同规模的样 
本，样本1的规模为，样本2的规模为 N 2 。 假定两个样本 
的协方差矩阵是一样的，既然 △ 与样本规模无关，那么两个 
样本中的△也是一样的。然而，由于样本规模不同，两个 
样本的卡方值也不相同。具体来讲 ， xf ：[)△， X 〖= 

( N 2 -1) A , 两者之间的关 系为： 

z 统计量与样本规模的关系也是类似的。 w 为任意一个 
具体的参数心抽样分布的方差为 [ A/(N — 1)]，△是与样本 
规模无关的值。同样，用同一个模型来估计两个不同规模的 
样本，样本1的规模为\，样本2的规模为 N 2 。 假定两个样 
本的协方差矩阵是一样的，既然△与样本规模无关，那么两个 
样本中的△也是一样的。然而，由于样本规模不同，抽样分布 
会不同。具体来说， o \ = A /( N ! — 1) , a \ = A /( N 2 — 1) ，两者 

之间的关系为： 



z 统计量的公式为 z = (£ / )/ cr ， oT 为假设中设定的常 
数，既然£与/都不受样本规模的影响•则样本规模对样 
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协方差结构横型 : USREL 导论 


本统计量的影 响为: 



其中， Zl 为样本规模为 N , 的样本统计量 ， Z * 2 为样本规模为 
N 2 的样本统计量。 

模型修正指标 

如果模型拟合较差，可以进行模型设定搜寻过程。例 
如，可以删除那些统计上不显著的参数，当然也可以加入某 
些参数。通过释放具有最大模型修正指标的参数，通常可以 
很好地改善模型的拟合度。在模型设定搜寻过程中，需要注 
意以下 几点: 第一，必须证明所考虑的模型都是可以识别的; 
第二，不能仅通过常规的 z 检验或卡方检验来选择模型，这 
是因为模型是由数据选择的，因此不能用同样的数据来检验 
模型; 第三，必须在理论和实际的指导下搜寻模型。即便是 
根据理论所设定的模型被拒绝了 •也不能为了增加模型的拟 
合度而删除或者增加某个参数。 

衍生参数 

计算由模型估计得岀的衍生参数通常是很有帮助的。 
可以通过估计的协方差矩阵计算潜变量的相关系数。例如， 

&和乞 之间的相关系数可以通过来计算， 
其中0为矩阵中的元素。为了评估潜变量的多个指标，也 
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可以计算其信度(关于信度的计算与解释详细参阅《验证性因 
子分析》)。对于模型中的结构部分，可以计算每个方程的确 
定系数。例如，对于预测 7 ，的方程，其确定系数为1^=(^一 
士 d / h (详细参考第3章中关于确定系数的讨论）。前文因 
子模型及结构方程模型中关于参数的讨论都可以直接应用 
到协方差结构模型的分析中。 

下面我们通过两个例子来说明上文的讨论。 

例 4.1: 参数估计与假设检验 

表 4.2 列出了关于心理失调的六个模型的 ML 估计结 
果。模型 Hi 与惠顿 （ Wheaton ，1978) 估计的模型最相 

近。 [1 G ] 虽然模型 M a ， M b ， 不太符合实际情况，但是它们 
对说明模型的设定和假设检验是很有帮助的。模型 M t .和 
M f 包含了误差项的序列相关，是对惠顿模型的扩展。 

模型 M a 假定所有方程中的误差项都是不相关的，即所 
有变量的误差都不相关，因此没有任何参数是相等的。模型 
M a 的卡方值为142.6,自由度为13•这表明该模型的拟合度 
很差。模型增加了两个合理的系数等同 限定： （1)/? 21 = 
仏 2 ，表明前期 SES 对后期 SES 影响的稳 定性； (2) A ^ 3 =入“， 
心理失调症状的因子负载在不同时期的尺度是一样的。增 
加了这两个限度后，自由度增加了 2个，卡方增加了 4.84。 
我们可以通过比较模型 M a 和模型 M b 的卡方差异检验这两 
个参数限定，即 Ho ：^ 21 =/? 42 , AJ 3 = A &。 卡方差异为 X 2 = 
X 2 b — X : =4.84, 自由度 df = df a — df b =2, 因此不能在95%的 
置信水平上拒绝原假设。 
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协方差结构 模型: LKREL 导论 


模型 H . 增加了参数心和仏 7 ，即关于同一变量的测量 
误差在不同时间点上是相关的。加入这两个参数的理由如 
下 :例如 ，时点2上影响心理失调症状（％)测量误差的因素 
很有可能也会影响时点3上心理失调的症状(％)的测量误 
差，即‘关0,同样，^也不等于0。释放这两个参数用掉了 
2个自由度，卡方差异减少了 16.5,因此可以拒绝4 =# 17 =0 
的原假设。 

模型 M a 和模型 Mi 、 分布用了相同的参数再现协方差结 
构，可是模型更好一些，因为它的卡方值比较小。可是由 
于 M a 和 M t . 并不是嵌套模型，所有无法检验 11.7 的卡方差 
异是否显著。可是基于实际情况的考虑，是一个更好的 
模型。 

模型 M H 在的基础上增加了同时点上误差相关的限 
定。释放参数少 23 和^ ^之后，模型减少了 2个自由度，卡方 
值减少了 41.2。与模型 M t . 相比，模型 M d 的拟合度有了显 
著的改善，可是整个模型的拟合度在统计上依然是不可接受 
的 ( X 2 =89.8, df = ll )， 这个拟合度与惠顿 ( Wheaton ，1978) 
所估计的模型是相似的。 

尽管我们是根据理论设定的模型 M d ，可是其拟合度并 
不好。 M d 中最大的修正指标是系数少 24 ，其修正指数为 
36.6。少 24 表明在时点2与时点3时预测 SES 方程的误差存 
在相关，这是合理的，因为模型中没有包含影响 SES 的变量。 
释放这一参数后，模型的自由度减少了一个，其卡方值减 
少了 44.5,显著地改善了模型拟合度。模型 M f 将这一限定 
扩展到不同时点时预测心理失调症状方程的误差项是相关 
的，尽管模型的拟合度并没有改善太多，但是在95%的置信 
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水平上是显著的。 

模型 M f 的自由度为 9, 卡方值为 40.8, 在统计上讲，其拟 
合度依然不能被接受。模型修正指标显示，如果增加参数 
仏 6 ，模型的卡方会改善 17.37。 然而，既然假定了 ^的协方差 
为0,释放^和之间的协方差是没有意义的。这说明当利 
用模型修正指标来进行参数设定时，必须注意如果所增加的 
参数不符合理论或实际，那么对模型拟合的改善也是没有意 
义的。 

协方差模型中的结构部分与结构方程模型是一致的（参 
见表 3.4 中的模型 M b )， 不同的是在模型 M d 中由时点2心 
理失调到时点 3 的 SES 的路径系数是显著的，这在表 3.4 的 
模型 M h 与惠顿的模型中都是不显著的。这条路径系数是很 
重要的，惠顿将其称作社会选择效应，即时点 2 和时点 3 之 
间的地位获得受到了心理失调的影响。这与社会因果效应 
是相反的，即由时点1的 SES 到时点 2 的心理失调的路径系 
数。此处的分析同时为社会因果效应(爲,）与社会选择效应 
(仏 3 )提供了证据，而惠顿只发现了社会因果效应。 

表 4.2 心理失调协方差结构的最大似然估计 

参数 模型 M a 模型 M b 模型 Me 模型 Md 模型 Me 模型 Mf 

yu 0.294 … 0.294 刪 0.294*** 0.294*** 0.294 … 0.294 糾 

/?2i 0.877 … 0.862*** 0.861 a ** # 0.860 a## * 0.894 a## * 0.893 a ** # 

y 2 i 0.026 0.031 * 0.031 * 0.031 * 0.021 0.021 

/?3i -0.010 刪 -0.010 -0.010*** 一 0.008*** -0.008 紛 -0.006*** 

y 31 - 0.002 - 0.002 - 0.002 0.000 — 0.000 — 0.000 

^42 0.844 *** 0.862 3 *** 0.861 8 *** 0.860 a ** # 0.894 a * # * 0.893 a ^ 

/? 43 0.052 0.160 -0.002 -1.495*** -0.846** -0.916* 

y 4 i 0.061*** 0.056*** 0.055*** 0.053*** 0.044** 0.045”* 

/?52 -0.000 0.000 -o.ooo -o.ooo -0.000 0.001 
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协方差结构 模型: LISREL 导论 


续表 


参数 

模型 M a 

模型 M h 

模型 Me 

模型 Md 

模型队 

模型 M f 

卢53 

0.601”* 

0.661 w 

0.619 … 

0.672 •“ 

0.649 … 

0.949 — 

/51 

- 0.002 

- 0.002 - 

- 0.002 

-0.000 

- 0.000 

0.000 

COR ( ^ . ^ 2 ) 

— 

一 

_ 

0.293 … 

0.269 *.• 

0.243*** 

COR ( h，（2 ) 

— 

一 

— 


-0.276 … 

—0.268 •“ 

COR ( C 5 • (2 ) 

— 

— 

一 

— 

一 

—0.389 … 

COR ( C ，⑺ 

— 

一 

— 

0.186 

0.181 

0.122 


0 . 211 *** 

0.244 b### 

0.268 b *** 

0.485 b "- 

0.443 b ** # 

0.503 b “* 

^75 

0.271 *** 

0.244 1 ，… 

0.268 bw 

0.485 b *** 

0.443 b * # * 

0.503 b *** 

REU173 » >»3) 

0.823 

0.746 

0.676 

0.367 

0.408 

0.361 

REL ( J73 . y \) 

0.250 

0.304 

0.331 

0.588 

0.546 

0.625 

REL ( 75 ， ) 

0.740 

0.792 

0.729 

0.413 

0.446 

0.392 

RELCays * y ?) 

0.407 

0.355 

0.394 

0.730 

0.657 

0.749 

C () R ( f 3 w 6 ) 

— 

— 

0.278 

0.465 … 

0.457 … 

0.469 … 

COR ( m “ 7 ) 

一 

— 

0.140 

-0.364 ” 

— 0.200 

-0.386 “ 

卡方值 

142.65 

147.49 

131.00 

89.84 

45.37 

40.79 

df 

13 

15 

13 

11 

10 

9 

P 值 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


注 ：罗马 字母表示该参数限定等同。 * 代表 p < 0 . 1 ，** 代表 p < 
0.05, *** 代表 p < 0.025。关于系数的解释可以参照图4.1。所有的系数 
都没有标准化，采用的是完全信息最大似然估计法。 


读者不应该将我们与惠顿的分析的差异看做对其的批 
评，本书在应用惠顿的模型时，做了许多简单的假定。这些 
简化是研究者为了设定模型所强加的，而它们过于简单。我 
们所发现的重大差异是以协方差模型的效力为代价的。由 
于这个模型的复杂性，在设定模型时，需要作出很多没有实 
际意义的假定，这些假定可能会直接影响到从模型中得出的 
实际结论。因为采用全信息模型，模型设定中的细微改变都 
有可能得出完全不同的结果。例如，从模型 H . 到模型 M d . 
社会选择效应(仏 3 )就经历了从不显著到显著的改变。同样， 
从模型 Ni 和模型，读者也可以看到改变其中一个参数 
会影响到模型中的其他参数。 
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模型中测量部分的重要意义在于保证了各自测量的信 
度。从模型 M a 到模型叭，心理失调的测量和％)是最 
稳定的。然而，在模型 M d 中，关于精神生理失调症状的测量 
是最稳定的(％和^ ) ，而且从模型到模型 M d ，测量误差 
的协方差也变得显著了。 

例 4.2: 参数估计与假设检验 

本例在例 3.3 的基础上为每个内生变量增加了两个测量 
指标（比较图 3.2 和图4.2)。在第3章中，我们讨论到由于影 
响研究对象期望的因素和影响其朋友期望的因素是对称的， 
因此其对应的参数也应该是相等的。例如，研究对象智力水 
平对其期望的影响 （ y l2 ) 应该与其朋友的智力水平对其朋友 
的期望的影响 （ y 25 ) 是相等的。同样，研究对象职业期望 (3^) 
在其期望(％)上的因子负载与其朋友的职业期望 （3^) 对其 
自身的期望( 7 2)的影响是相等的，即 = A ( 2 。 此设定的结 
果在表 4.3 的模型 M a 中，该模型的自由度为24,卡方值为 
30.6,表明模型的拟合度是很好的。 


表 4.3 同辈对期望影响的协方差结构模型 


参数 

模型 M a 

模型 M b 

/?12 

0 .202*"* 

0 .201 **** 

711 

0.166 b … 

0.170 b — 

Y\2 

0.307”“ 

0.310 c *.* 

713 

0.233 d … 

0.235 dw 

714 

0.071 … 

0.074 一 

p2\ 

0 .202° 

0 .201”“ 

723 

0.071”“ 

0.074 … 

/24 

0.233 d .** 

0.235— 
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续表 


参数 

模型 M a 

模型 Mb 

725 

0.307 c *** 

0.310 c ** # 

/26 

0.166 b *“ 

0.170 b “* 

COR(a ， 0 ) 

-0.065 

-0.162 

x y 

a 21 

0.936 f *，* 

0.908 f '** 

A 4 y 2 

0.936 f *** 

0.908 f * # * 

REL ( , yi ) 

0.450 

0.475 

REL(^i ， y2 ) 

0.394 

0.392 

REL( 72 • y3) 

0.464 

0.490 

REL( 72 , y4) 

0.406 

0.403 

CORG 2 , 

— 

0.267*** 

COR( n » rs ) 

— 

0.065 

卡方值 

30.63 

16.10 

df 

24 

22 

P 值 

0.17 

0.81 


注 :罗马 字母表示该参数限定等同。 * 代表 p <0.1, ** 代表 p < 
0.05, *** 代表 p < 0.025。 关于系数的解释可以参照图 4.2, 所有的系数 
都没有标准化，采用的是完全信息最大似然估计法。 


该模型结构部分的结果与例 3.3 中结构方程 t 旲型的结果 
是一样的。在例 3.3 中由观测变量导致的变化，在这里就变成 
了由通过两个指标所反映的潜变量所导致的变化。对于研究对 
象和其朋友而言，都是教育期望的指标比职业期望更为可靠。 

模型同时测量了研究对象与其朋友的教育期望，考虑到 
同一个指标用了两次，很有可能两个变量的误差是相关的, 
即心參0。 同样，职业期望的测量误差也是相关的， 

0。模型 M b 添加了这两个参数，增加了这两个参数后，模型 
M b 的自由度减少了 2个，卡方减少了 14.5,表明模型拟合有 
了很大的提高。与模型 M a 相比，模型 M h 中的其他参数并 
没有太大的改变，表 4.3 显示模型 M b 的自由度为22,卡方值 
为 16.1。// 


第 4 章协方差结构模型 


第8节 | 小结 



协方差结构模型是很复杂的•本书的介绍比较简单，因为 
我们没有介绍关于其数学模型的发展。掌握这个模型的最简 
单的方法就是学习并且应用。下面我们为读者推荐一些关于 
协方差结构模型的书本或文章 : Alwin Jackson , 1982； Ba - 
gozzi , 1981 ； Bielby et al . ， 1977 ； Bynner , 1981; Dalton ， 
1982； Joreskog , 1974， 1978； Jundd Milburn ，1980； Ken ¬ 
ny , 1979； Kessler Greenberg ，1981 ； Krehbiel Niemi , 
1982； Long , 1981 ； Sullivan Feldman ， 1979 。 LISREL 4.0 
和 LISREL 5.0 的手册也包含了大量有用的应用 ( J 6 resk 0 g & 
Sorbom , 1978， 1981) 0 
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因子模型、回归与相关模型、多元指标模型、二阶因子 
模型等都可以被整合到协方差结构模型分析的框架下。可 
是，尽管协方差结构模型有很强大的效力和灵活性，我们需 
要注意这个模型也存在着一些局限：第一，到目前为止，还 
没有关于跨组的模型比较。第二，尽管在实际研究中有许 
多变量是分类或定序的， ML 和 GLS 估计方法仍假定观测 
变量服从正态分布。第三，模型假定线性关系，对参数的限 
定要么为0.要么相等。有大量的研究者致力于扩展这方面 
的研究，例如，本特勒、约瑞斯科、麦克唐纳 （ McDonald )、 穆森 
( MMhen )、 索尔波姆、威克斯等。《测量心理学》、《多元行为 
研究》、《英国统计与数学心理学期刊》、《计量经济学》等刊物 
会定期发表该模型研究的最新进展。关于协方差结构模型 
最新的研究进展的详细讨论超越了本书的内容，下面我们为 
读者提供协方差结构模型最新发展的简单介绍。 

1. 分组比较。检验组间差异是很多研究的基本问题，比 
如变量: r 的影响对男性和女性的影响是否相同？控制组与实 
验组之间是否有差异？约瑞斯科和索尔波姆在 LISREL 5.0 和 
6.0 中整合了协方差结构模型分析的组间比较程序。此外，本 
特勒、威克斯和麦克唐纳的讨论也包含了组间比较。组间差 
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异比较要求同时对每个组估计其协方差结构，并且在组间比 
较各种参数限定。关于组间差异比较的应用可以参考索尔 
波姆和马吉德森的论著 ( Sbrbom ， 1982； Magidson ， 1977)。 

2. 离散数据分析。 ML 和 GLS 的估计方法假定观测变量 
服从正态分布。如果数据是非正态的，建议用 ULS 估计法。 
如果变量是分类的或者定序的，可以假定它们反映了一个连 
续的不可观测的变量，在用 ULS 估计法时，输入 Polychoric 、 
Tetrachoric 和 Polyserial 相关系数矩阵（参考 Joreskog 8^ 
Sorbom , 1981:第6章）。约瑞斯科和索尔波姆展示了如何 
通过 LISREL 5.0 来分析离散数据。 

3. 更为复杂的结构与限制。我们前文曾指岀协方差结 
构模型可以通过参数的任何函数来解释观测变量的协方差 
(假定模型是可以识别的）。麦克唐纳 （ McDonald , 1978) 和 
约瑞斯科 ( hreskog ，1978) 讨论了各种函数形式，尽管这些 
函数形式很难应用到实际分析中。本特勒和麦克唐纳发展 
了一些不是那么普遍但是具有重要应用价值的模型。这些 
模型的优势在于包含更为广泛的参数限定形式及更为复杂 
的测量模型。 

4. 协方差结构模型的软件。无论是基本的协方差结构 
模型还是较为复杂的协方差结构模型，都需要比较复杂的软 
件来估计。确实，协方差结构模型的每次重大突破都伴随着 
统计软件上的进步，因此有必要简单介绍一些可以用来估计 
协方差结构模型的软件程序。 

第一个用来估计协方差结构模型的软件是由约瑞斯科 
和范西洛 (Rreskog & van Thillo , 1972) 在他们进行教育测 
验时写作的 LISREL 。 这个程序的第一个版本含有冗余的控 
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制程序而且不能估计标准差•其优势在于它是第一个可以估 
计协方差结构模型的软件而且是免费的 。 LISREL 2.0 包含 
了标准误的计算，读者需要注意，这个版本里有些严重的错 
误。最新的版本是由国际教育服务发布的 LISREL 5.0®，这 
个版本简化了许多模型设定、 ULS 和 ML 估计、标准误的估 
计、模型修正指标、初始值设定等的语言命令，而且还可以进 
行组间比较及定序变量的分析。 

MIIS 是由国家精神健康机构的罗纳德 • 舍恩伯格 
(Ronald Schoenberg ) 用 R 编写的，而且是免费的。 MILS 可 
计算标准误、 GLS 和 ML 估计、组间比较，而且扩展了一些包 
含乘法模型的协方差结构模型。 MILS 的缺点在于其命令操 
作比较难•不提供初始值计算，而且没有提供一些有用的统 
计指标（比如模型修正指标）。其优势在于可以评估参数解 
的稳定性及计算间接效应的标准误。 

由彼得 • 本特勒 (Peter Bender ) 发展的 EQS ， 包含了许多 
LISREL 5.0 的特征，如可以进行组间比较，计算初始值，提供 
ML 好 ULS 的估计。此外， EQS 还可以允许复杂的参数限 
定，计算 GLS 估计，以及估计由本特勒和威克斯 （Bentler & 
Weeks , 1979) 发展出来的更为普遍的协方差结构模型。 

在某些方面，由麦克唐纳发展的 COSAN 是目前最为普 
遍的程序，它允许研究者写作 FORTRAN 子程序来估计非 
标准化的模型。 

考虑到协方差结构模型及估计软件的复杂性，研究者极 
有可能在写命令时出错，这不仅会导致错误的估计，而且也 


①目前最新的版本是 LISREL 8.7。 译者注 
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极为耗时。经验表明，正确的模型通常比错误的模型更为省 
时。除非研究者非常熟悉估计软件，重新估计已经发布的文 
章对于理解估计程序是很重要的。本书中的例子所用到的 
数据都包含在附录里•读者可以运用这些数据重新估计本书 
中的范例。 
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范例中的协方差/相关矩阵 


表 A 惠顿模型的相关矩阵 (N=603) 


变量 

FSES 

SES1 

SES2 

SES3 

PSY2 PHY2 PSY3 

PHY3 

FSES 

1.000 








SES1 

0.257 

1.000 







SES2 

0.248 

0.882 

1.000 






SES3 

0.264 

0.827 

0.863 

1.000 





PSY2 

-0.072 

— 0.166 

一 0.096 

— 0.089 

1.000 




PHY2 

一 0.038 

— 0.104 

— 0.016 

—0.030 

0.454 

1.000 



PSY3 

-0.092 

— 0.120 

一 0.088 

— 0.006 

0.526 

0.247 

1.000 


PHY3 

— 0.013 

-0.139 

-0.011 

— 0.005 

0.377 

0.309 

0.549 

1.000 

S.D. 

19.98 

22.82 

22.85 

22.69 

1.45 

0.555 

1.38 

0.503 


注 :变量 的识別(括号中的标注表明分別对应例 3.2 和例 4.1)： 

FSES 为父亲的社会经济地位(心， ：n); 

SES 1 为研究对象时点1时的社会经济地位( 71 ， yx ); 

SES 2 为研究对象时点2时的社会经济地位( 72 , ： y 2 ); 

SES 3 为研究对象时点3时的社会经济地位( 7< , ^); 

PSY 2 为时点2时心理失调症状的数目（ 73 , ^ 3 ); 

PSY 3 为时点3时心理失调症状的数目 （ none ， yO ； 

PSY 2 为时点2时精神生理失调症状的数目 （ V5 , ys)i 
PSY 3 为时点3时精神生理失调症状的数0 ( none . : y 7 )。 

资料 来源： Blair Wheaton. ** Sociogenesis of Psychological Disorder* w 
American Sociological Revieu^ ♦ vol.43，pp.383 — 403( 表 2)。 


。(1 懈 ) ZCOI—elldd *5S *S070 

c < -so^£^ax< cisacc^ ;^S-«X.CT< -csso.c.o 

* ( f<v ) 挪萚 s^c^^s 鉍友炽 fe^Pusvd-J 

---H 

K 一 -sr@lste^s 礙友饫 fe^dsv3->i 
乂 zk ： *auou)© 萚， 4p”f 运礙灰保宿武 dsv(s 
K? *<Dsm 纥颂彻龙运礙友饫窖武 sasi 
乂 二：一 u/)©l 萚 S 电 CX 礙 fe 饫串武 dsvcli 
乂？ *<D)aI S 槪友炽宙积 aTH 
"(rH 运运 ecof 在但友菡汆菸蝌坦遊 §*&- 妒靼 ) 系琪§_$.:迕 


F-EASP 










1.000 0 

F-( )ASP 









1.000 0 

o 

Z 

• 

O 

F-SES 








1.000 0 

0.360 7 

0.410 5 

F-PASP 







1.000 0 

-0.043 8 

0.198 8 

0.278 4 

g 

u- 






1.000 0 

0.208 7 

0.295 0 

0.500 7 

0.519 1 

R-EASP 





1.000 0 

0.286 3 

0.070 2 

0.240 7 

0.327 5 

0.366 9 

R-OASP 




1.000 0 

0.624 7 

0.299 5 

0.076 0 

0.393 0 

0.421 6 

0.326 9 

R-SES 



1.000 0 

0.324 0 

0.404 7 

0.2302 

0.093 1 

0.270 7 

0.278 6 

0.305 4 

R-PASP 


1.000 0 

0.048 9 

0.213 7 

0.274 2 

0.078 2 

0J147 

0.018 6 

0.083 9 

0.1124 

a 

1.000 0 

0.183 9 

0.222 0 

0.410 5 

0.404 3 

0.335 5 

0.102 1 

0.186 1 

0.259 8 

0.290 3 

B 

R-IQ 

R-PASP 

R-SES 

R-OASP 

R-EASP 

F-IQ 

F-PASP 

F-SES 

F-OASP 

F-EASP 
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注释 


[1] //表示范例已经结束。 

[2] 更精确地说，这个条件意味着，通过重新排序方程 ,B 可以简化为三角 
矩阵。 

[3] (rXc) 维的矩阵 X 的秩就是 X 子矩阵中最大的非奇异矩阵.例如，矩阵 



•1 2~ 

的秩为2,因此 X 是可 逆的； 再例如， X = 

12 1' 

X = 

Q A 

3 4 3 




1 1 1 


秩为2,因此X是不可逆的，但是 



是矩阵 X 的最大子矩阵而且 


是可逆的。 

[4] 这里的一般最小二乘法是指单个方程的估计方法.有时特指 Airken 的 
最小平方。 

[5] 这种估计方法与计量经济学是不同的。 ULS 和 GLS 的估计值与文献 
中常见的是不一样的。 ML 对应完全信息最大似然，这与马林沃德 
(Malinvaud, 19 7 0> 的估计是相似的。 

[6] 恰好识别模型的自由度为0,但并不是自由度为0的都是恰好识别模 
型，不能识别的模型的自由度也可以为0。 

[7] 关于因果模型中的直接效应、间接效应和总效应的详细处理可以参考 
Wonnacott 8 •- Wonnacott ， 1981:194—207； Alwin Hauser, 1975； 
Duncan ， 1975。 

[8] 惠顿指出，尽管在追踪数据模型中，通常会施加等同限定，但是在这里 
并不适合。例如亨内平 （Hennepin) 地区在时点2和时点3之间经历了 
工业化，然而在时点〗到时点2之间并没有工业化。因此，社会经济地 
位 (SES) 的稳定性可能会改变。然而，为了教学，我们采用这种等同限 
定是很有帮助的。 

[9] 这些结果并不仅适合协方差结构模型，它们也可以应用到统计推断的 
其他问题中。 

[10] 我们的设定与惠顿的设定有两个主要差异 :第一 ，惠顿假定关于 SES 的 
测量是不完美的•我们假定其测量是完美的；第二，我们假定不同时间 
上早期社会经济地位对后来社会经济地位影响的稳定性.而惠顿并没 
有采用这一假定(注释[8])。此处的模型修正是出于教育目的.而非实 
际情况。我们建议读者阅读原文。 




参考文献 


1() ( J 


参考文献 


ALWIN, D.F. and R.ML HAUSER (1 975 V*The decomposition of effects in 
path analysis. American Sociological Review 40 ： 37 — 47. 

ALWIN, D.F. and D.J. JACKSON(1982) “ 丁 he statistical analysis of Kohn's 
measures of parental valuevS,” pp. 197 — 223 in H. Wold and K. Joreskog 
(eds.) Systems Under Indirect Observation• New York ： Elsevier North- 
Holland. 

BAGOZZI, R.P. (1981) “An examination of the validity of two models of 
attitude.” Multivariate Behavioral Research 16 ： 323 — 359. 

BENTLER, P.M. (1980) “Multivariate analysis with latent variables ： causal 
modeling.” Annual RevieiL' of Psychology 31:419 — 456. 

BENTLER, Peter M. (1989) EQS : Structural Equations Program Manual. 
Ix>s AngelevS ： BMDP Software. 

- (1976) “Multistructure statistical model applied to factor analysis.” 

Multivariate Behavioral Research 11:3 —— 25. 

- D. G. BONETT ( 1980) “Significance tests and goodness-of-fit in the 

analysis of covariance structures.” Psychological BuLietw 88 ： 588 — 606. 

BENTLER, P.M. and D.G. WEEKS(1979) “Interrelations among models 
for the analysis of moment structures.” Multivariate Behavioral Re¬ 
search 14 ： 169 — 185. 

BIELBY, W.T., R. M. HAUSER and D.L. FEATHERMAN(]977) "Re¬ 
sponse errors of nonblack males in models of the stratification process.” 
Journal of the American Statistical Association 72:723 —— 735. 

BLAU ， P. and O. D. DUNCAN ( 1967) American OccupationaL Structure. 
New York ： John Wiley. 

BOCK. R.D. and R.E. BARGMANN(1966) “Analysis of covariance struc- 
tures.” Psychometrika 31:507 — 534. 

BOLLEN，Kenneth A. (1989) Structural Equations with iMtent Variables. 
New York ： John Wiley. 

BOOMSMA, A. (1982) “The robustness of LISREL against small sample si¬ 
zes in factor analysis models，’’ pp. 149 — 173 in H.Wold and K.Joreskog 
(eds.) Systems Under Indirect Observation. New York ： Elsevier North- 
Holland. 

BROWNE, M. W. ( 1974 ) Generalized least-squares estimators in the 



11() I 协方差结构 横型: USREL 导论 


analysis of covariance structures.” South African Statistical Journal 8 ： 
1—24. 

BYNNER, J.(1981) “Use of LISREL in the solution to a highei^orcler factor 
problem in a study of adolescent self-images. M Quality and Quantity 
15j 523—540. 

DALTON, R.J.(1982) “The pathways of parental socialization/’ American 
Politics Quarterly 10 ： 139 — 157. 

DUNCAN» O. D. ( 1975) hit rod uct ion to Structural Equation Models. New 
York ： Academic. 

- A. C). HALLER, and A. PORTES ( 1971 ) Peer influences on 

aspirations ： a reinterpretation. M pp.219 -241 in H.M. Blalock. Jr.(ed) 
Causa! Models in the Social Sciences. Chicago ： Aldine. 

FISHER, F. M. ( 1966) The IJev/ificarion Problem in Econometrics. New 
York ： McGraw-Hill. 

F()X ， J.C1980) “Effect analysis in structural equation models. Sociological 
Methods ami Research 9 ： 3 — 28. 

GOLDBERGER» A.S.( 1971) **Econometrics and psychometrics ： a survey of 
communalities •” Psychometrika 36:83 — 107. 

- O.D. DUNCAN (1973) Structural Equation Models in the Social Sci- 

ertces. New York ： Seminar. 

GRAFF. J. and P.SCHM1DT( 1982) **A general model for decomposition of 
effects. M pp. 197 — 223 in H.Wold and K.Joreskog(eds.) Systems Under 
Indirect Ofjservation. New York ： Elsevier North-Holland 

GUILFORD, J. P. and Ii FRUCHTER (1978) Fu?iJamental Statistics in 
Psychology and Education. New York ： McGraw-Hill. 

HANUSHEK ， E.A. and J.E. JACKS()N( 1977) Statistical Methods for So¬ 
cial Scientists. New York ： Academic. 

JORESKCXj. K.G.(1978) 44 Statistical analysis of covariance and correlation 
matrices.” Psychometrika 43:443 — 477. 

- (1974)** Analyzing Psychological data by structural analysis of covariance 

matrices,” pp. 1 —— 54 in R. C. Atkinson et al. (eds.) Contemporary 
Developments in Mathematical Psychology * vol. 2. San Francisco ： 
Freeman. 

- (1973) “A general method for estimating a linear structural equation 

system，’’ pp.85 — 112 in A.S. Goldberger and O.D. Duncan (eds.) Struc¬ 
tural Equation Models in the Social Sciences. New York ： Seminar. 



参考文献 


111 


- A. S. GOLDBERGER ( 1972 ) Factor analysis by generalized least 

squares.” Psychometrika 37:243 — 260. 

K)RESKOG，K.G. and D.S()RBOMC1981) USREL V.User's Guide. Chic- 
gao： National Educational Resources. 

- (1978) LISREL IV. User's Guide. Chicago： National Educational Re¬ 
sources. 

- (1976) LISREL III : Estimation of Linear Structural Equation Sys¬ 
tems by Maximum Likelihood Methods. User's Guide. Chicago： Inter¬ 
national Educational Services. 

J()RESKOG, K.G. and Mvan THILLO(1973) LISREL： A General Com¬ 
puter Program for Estimating a Linear Structural Equation System In¬ 
volving Multiple Indicators of Unmeasured Variables. Research report 
73 — 5. Department of Statistics，Uppsala University, Uppsala Sweden. 

J()RESKOG, Karl G., and D.S()RB()M( 1988) LISREL 7： A Guide to the 
Program and Applications. Chicago： SPSS, Inc. 

- Mvan Thillo( 1972) LISREL : A General Computer Program jur Ksti- 

mating a Linear Structural Equation System Imvlvin^ Multiple Indicators of 
Unmeasured Variabies. Princeton, NJ： Educational Testing Service. 

JUDD, C. M and M. A, MILBURN (1980) “The structure of attitude 
systems in the general public.” American Sociological Review 45： 
627—643. 

KENNY♦ D. A.(1979) Correlation and Causality. New York： John Wiley. 

KEESLING, W. ( 1972) M Maximum likelihood approaches to causal flow 
analysis.” Ph.D. dissertation University of Chicago. 

KESSLER，R. C. and D. F. GREENBERG (1981) Linear Panel Analysis. 
New York： Academic. 

KMENTA, J.( 1971) Elements of Econometrics. New York： Macmillan. 

KREHBIEL* K. and R.G. NIEMK1982) 44 A new specification and test of 
the structuring principle. v Paper delivered at the annual meeting of the 
American Political Science Association. 

LEE. S.Y.C1977) “Some algorithms for covariance structure analysis.” Ph. 
D. dissertation* University of California* Los Angeles. 

LONG. J.S.( 1981) “Estimation and hypothesis testing in linear models con¬ 
taining measurement error,” pp.209 — 256 in P.V. MarsdenCed.) Linear 
Models in Social Research. Beverly Hills，CA： Sage. 

McDONALD，R.P.(1980) “A simple comprehensive model for the analysis 


112 


协方差结构 模型： USREL 导论 


of covariance structures： some remarks on applications.” British 
Journal of Mathematical and Statistical Psychology 33： 161 — 183. 

McDONALD，R.P.(1978) “A simple comprehensive model for the analysis 
of covariance structures.” British Journal of Mathematical and Statis¬ 
tical Psychology 31:59 — 72. 

MAGIDSON，J.(1977) “Toward a causal model approach for adjusting for 
preexisting differences in the nonequivalent control group situation： 
general alternative to ANCOVA. ” Evaluation Quarterly 1:399 — 419. 

MALINVAUD, E. (1970) Statistical Methods of Econometrics. New York： 
Elsevier North-Holland. 

SCHOENBERG* R. (1982) MILS : A Computer Program to Estimate the 
Parameters of Multiple Indicator Linear Structural Models. Bethesda, 
MD： National Institutes of Health. 

SORBOM, D.(1982) * 4 Structural equation models with structured means，” 
pp. 183 — 195 in H. Wold and K. Joreskog (eds.) Systems under Indirect 
(y)servation. New York： Elsevier North-HollanA 

- (1975) “Detection of correlated errors in longitudinal data. M British 

Journal of Mathematical and Statistical Psychology 28:138 — 151. 

SULLIVAN♦ J. L and S. FELDMAN (1979) Multiple Indicators. Beverly 
Hills, CA： Sage. 

THEIL, H. (1971) Principles of Econometrics. New York : John Wiley. 

WHEATON. 11(1978) “The sociogenesis of psychological disorder.” Amer¬ 
ican Sociological Review 43： 383 — 403. 

—— B.MUTHEN，D.ALWIN, and G.SUMMERS(1977) a Assessing relia¬ 
bility and stability in panel models, M pp.84 — 136 in D.R HeiseCed.) So¬ 
ciological Methodology . 1977. San Francisco: Jossey-Bass. 

WILEY* D. E. (1973) “The identification problem for structural equation 
models with unmeasured variables,” pp.69 — 83 in A.S. Goldberger and 
O.D. Duncan(eds.) Structural Equation Models i?t the Social Sciences. 
New York： Seminar. 

WONNACOTT，R.J. and T. H. WONNACOTTC1979) Econometrics. New 
York： John Wiley. 

WONNACOTT，R.H. and RJ. WONNACOTTC 1981) Regression : A Sec¬ 
ond Course in Statistics. New York： John Wiley. 



译名对照表 


113 


译名对照表 


Confirmatory Factor Analysis ( CFA ) 验证性因子分析 

Covariance Structure Models ( C^SM ) 协方差结构模别 


endogenous variables 
estimation 
exogenous variables 
factor loading 
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full information techniques 
Gereralized Least Squares ( GLS ) 
goodness of fit 
identification 
latent variables 

limited information techniques 

Maximum LikehoodCML ) 

measurement model 

observed variables 

order condition 

panel mode ) 

rank condition 

recursive model 

simultaneous equation systems 

specification 

structure equation model 
Unweighted Least Squares ( ULS ) 


内生变量 
估计 

外生变量 
因子负载 
完全信息技术 
一般最小二乘法 
拟合度 
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识别 
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潜变量 

有限信息技术 
最大似然估计 
测量摸型 
观测变量 


阶条件 








追踪数据模型 
秩条件 
递归模型 
联立方程系统 
设定 

结构方程模型 
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